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Summary 

A new procedure for analyzing multicomponent 
positron lifetime spectra in polymers has been devel- 
oped. It requires initial estimates of the lifetimes and 
intensities of various components, which are read- 
ily obtainable by a standard spectrum stripping pro- 
cess. These initial estimates, after convolution with 
the timing-system-resolution function, are then used 
as the inputs for a nonlinear least-squares analysis 
to compute the estimates that conform to a global- 
error minimization criterion. The convolution inte- 
gral uses the full experimental resolution function, 
in contrast to the previous studies in which analyt- 
ical approximations of it were utilized. These con- 
cepts have been incorporated into a generalized com- 
puter program for analyzing positron lifetime spectra 
(PAPLS) in polymers. The validity of this program 
has been tested by using several artificially generated 
data sets. These data sets were also analyzed with 
the widely used POSITRONFIT program. In almost 
all cases, the PAPLS program gives closer fits to the 
input values. The new procedure has been applied 
to the analysis of several lifetime spectra measured in 
metal-ion containing Epon 828 samples. The results 
are described in this report. 

Introduction 

Positrons are slowed down quickly (refs. 1 to 4) 
upon entering a condensed medium such as poly- 
mers. After thermalization, they may annihilate as 
free positrons or after being trapped in microvoids 
or at defect sites on molecular chains. Some of the 
trapped positrons also loosely bind with the molec- 
ular electrons to form positronium (Ps) atoms. (See 
refs. 3 to 8.) Depending on the final states of the 
positrons they annihilate at different times. The life- 
time of the positron can thus provide valuable infor- 
mation about the electronic and physical structure of 
the bulk material and the various defects or impuri- 
ties in it. 

Positron lifetime studies make use of radioactive 
sources that emit a gamma ray almost simultane- 
ously with the positron. The widely used Na 22 source 
emits a 1.28-MeV photon within 3 psec of the emis- 
sion of the preceding positron. The detection of this 
gamma photon marks the zero time (reference time) 
for a lifetime measurement, since the thermalization 
time for the positrons in polymers is much less than 
their subsequent lifetime and can usually be ignored. 
(See refs. 9 to 11.) The annihilation of the positron 
is signalled by the detection of one of the two anni- 
hilation photons produced in the [e + + e _ — + 2 hv 
(511 keV each)] process. The time interval between 
these two signals is a measure of the positron lifetime 


aad is commonly measured by fast timing systems. 
The two widely used timing systems are described in 
references 12 to 15. 

A typical lifetime spectrum of positrons annihi- 
lating in polymers (refs. 7 and 16 to 19) is expected 
tu consist of at least three exponentially decaying 
components, each characterized by a mean lifetime 
a; id a relative probability amplitude. The three com- 
ponents are short-lifetime (resulting from prompt 
positron and parapositronium decays), intermediate- 
li 'etime (resulting from trapped positron decays), 
aad long-lifetime (resulting from orthopositronium 
d 'cays). Quite often, the lifetime spectra in polymers 
have been analyzed by considering only two compo- 
nents, representing the mean lives of the “shorter” 
and “longer” components, for reasons of simplicity 
and/ or expediency. (See refs. 7 and 20 to 28.) How- 
e' er, a detailed analysis that provides at least three 
c< mponents is preferable, since the annihilations of 
all the positron groups are dependent on their lo- 
Cc 1 atomic environments and can thus provide more 
detailed information about the host-material defect 
st ucture. 

Observed spectra are subject to a finite time res- 
olution of the lifetime measurement system. Data 
events which would theoretically fall into a certain 
time interval (channel) actually exhibit a near Gauss- 
ia i distribution over many channels depending on 
the timing-system resolution. In previous studies, 
the timing-system- resolution function has been ap- 
proximated by an analytical expression for reasons 
of simplicity and convenience of analysis. However, 
the resolution function does not lend itself to a com- 
plete description by a reasonable analytical expres- 
sion. It would therefore be preferable to use the full 
ex oerimental resolution function for deconvolving the 
lifetime spectra. 

The purpose of this report is to develop a pro- 
cedure for resolving a multicomponent experimental 
sp *ctrum into individual lifetime components while 
imposing as few constraints and assumptions as pos- 
sible. This procedure has been incorporated into 
a i omputer program for analyzing positron lifetime 
sp *ctra (PAPLS) in polymers. This program has 
be >n used to analyze lifetime spectra observed in 
m< tal-ion containing Shell Epon 828 epoxy samples, 
an 1 the results are compared with those obtained 
wi h the previously developed analysis techniques. 
This procedure is adaptable for other types of con- 
densed matter within which positronium is likely to 
be formed. 


Problem Statement 

The lifetime spectra of positrons annihilating in 
a molecular solid consist of several components. It 
is often convenient to group the several lifetimes into 
three broad categories. (See ref. 29.) The short- 
est lifetime component (0.1 to 0.5 nsec) includes 
prompt free-positron decay and parapositronium de- 
cay. The intermediate lifetime component (0.5 to 
0.9 nsec) is believed to arise from the annihilation 
of trapped positrons. The longest lifetime compo- 
nent (1 to 4 nsec) arises from pick-off annihilation 
of orthopositronium atoms. Each of the annihila- 
tion processes contributes a decaying exponential, 
with a characteristic decay constant, to the spec- 
trum. In the previous studies of positron annihila- 
tion spectroscopy (PAS) in molecular solids (refs. 2, 
7, and 21 to 28), it has often been convenient to 
resolve the experimental lifetime spectra into two 
components. This procedure has been justified on 
the basis of the arguments that the long-life com- 
ponent, which results from the pick-off annihilation 
of the triplet positronium and is the most sensitive 
indicator of the host-material properties, can be eas- 
ily deconvolved from the rest of the spectrum; the 
short and intermediate lifetime components are then 
lumped together into a compound short-life compo- 
nent. It would, however, be more appropriate if all 
the positron groups could be deconvolved (refs. 16 
to 19) particularly since the other groups are also 
dependent on the atomic environment of the trapped 
positrons and can provide some very useful informa- 
tion about the host defect structure. 

For a three-component system, the lifetime spec- 
trum can be written as the sum of three exponentials 
as follows: 

3 

n (0 = ^i e + B 

i= 1 

= A ie ~ Xlt + A 2 e~ x * t + A 3 e~ X3t + B (1) 

where and A t represent zero-time amplitudes and 
decay constants, respectively, and B is a constant 
background. 

When experimental (n(£) versus t) data are dis- 
played on a semilog format, the spectrum appears to 
be a sum of three straight lines — each representing 
an exponential decay. The three groups are easily 
discernible if Aj, A2, and A3 are sufficiently differ- 
ent from each other and if their relative probabilities 
are comparable. In that case, the spectrum can be 
decomposed into individual components rather eas- 
ily. Sometimes, though, the A^ values are not suf- 
ficiently different, and the longer lived components 


are weaker; then, the spectrum analysis is a difficult 
problem. 

The theoretical spectrum defined by equation (1) 
must be modified for the finite time resolution of the 
lifetime measurement system. Data events which 
would theoretically fall in a certain time interval 
(channel) actually exhibit a near Gaussian distri- 
bution whose effective full width at half maximum 
(FWHM) equals the timing-system resolution. The 
form of the time-resolution function can be deter- 
mined by measuring a prompt spectrum, such as that 
obtained from detecting the near-simultaneous 1.17- 
MeV and 1.33-MeV gamma photons from a Co 60 
source, with the same experimental settings for elec- 
tronics that are used in the actual lifetime experi- 
ment. Zero time is determined by the centroid of 
this prompt spectrum. Incorporation of the finite- 
resolution function into equation (1) modifies the life- 
time spectrum as follows: 

/ +00 

R{t - t')n{t') dt' (2) 

-OO 

where R(t — t f ) represents the timing-system 
resolution-function and t f is the time variable, which 
differs from t by an integral number of time channels. 

In previous studies, the resolution function has 
been assumed to be Gaussian; therefore, a direct an- 
alytical deconvolution of the observed data was pos- 
sible. This approximation, however, tended to in- 
troduce systematic errors and to give poorer fits to 
the lifetime values. To improve on this situation, 
several authors have used modified Gaussian repre- 
sentations of the resolution function. Lichtenberger 
et al. (ref. 30) used a double-sided exponential func- 
tion, whereas Hall et al. (ref. 31) used a side-sloped 
function to represent the resolution function. Eldrup 
et al. (refs. 7 and 8) and Kirkegaard et al. (ref. 32), 
on the other hand, assume the time-resolution func- 
tion to be a sum of three Gaussian functions rather 
than a single Gaussian. Use of these various forms 
of analytical approximations improves the quality of 
the fits without complicating the solution of the con- 
volution integral. 

However, it has long been recognized that the 
experimental resolution functions of the fast timing 
systems used in PAS studies have slowly decaying 
characteristics which do not admit to simple ana- 
lytical representations. To overcome the limitations 
imposed by imperfect analytical representation of 
the resolution function, actual experimental resolu- 
tion functions are used as the smearing functions. 
An additional argument in support of this procedure 
arises from the fact that each timing system has its 
own unique resolution function. Even though this 
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approach makes the analysis procedure slightly more 
tedious, it is better suited for analyzing experimen- 
tal spectra with inadequate statistics, as is often the 
case with the more promising types of polymers such 
as PMR-15, LARC-TPI, and other high-temperature 
polyimides. 

The use of the actual experimental resolution 
function, rather than an approximate analytical 
representation of it, constitutes the major point of 
difference between the present and previous investi- 
gations. However, because of the numerical repre- 
sentation of the resolution function, it is not possible 
to analytically convolve the instrumental resolution 
effects into the theoretical spectrum. They have to 
be incorporated numerically as illustrated in the fol- 
lowing example. 

Assume that a lifetime spectrum has the following 
three components: 

r\ = 457 psec; h = 79.7% 

72 = 917 psec; I 2 = 5.3% 

73 = 2300 psec; 1 3 = 15.0% 

where r x is the lifetime of the ith component and I x is 
the intensity of the ith component. Also assume that 
the timing-system resolution is 250 psec, the time 
channel width is 61.23 psec, and the total counting 
time is such that the peak channel count is 10 5 and 
the background count is 10. Figure 1 shows the 
histogram that is observed with a perfect system. 
Random statistical fluctuations have been included 
in this calculated histogram. A timing system with 
a finite resolution will smear this histogram into 
a spectrum shown in figure 2. The total counts 
which fall in any one channel for a perfect system 
are spread over an experimental distribution defined 
by the Co 60 spectrum. This smearing process is 
repeated for each time channel in the theoretical 
histogram. Finally, the sums of all redistribution 
counts in the individual channels are recorded to 
give the computed spectrum (fig. 2) that would be 
observed with a finite-resolution timing system. 

Close comparison of figures 1 and 2 shows that 
the first 10 or so channels are severely disturbed by 
the finite-resolution effects. Channels higher than 
10 are essentially undisturbed, which indicates that 
the counts gained from the neighboring channels 
neutralize the losses to the neighboring channels that 
resulted from the finite-resolution effects. Figure 2 
also shows that the peak in the synthesized smeared 
spectrum is about 100 psec to the right of the time 
zero ( t z ) (reference) channel. This information can 
also be used to help locate the zero of the time 
scale in subsequent analyses. Figure 3 shows an 
actual experimental lifetime spectrum of positrons 


in an epoxy target. Obviously, the spectra shown in 
figures 2 and 3 are similar. 

In the following sections, a procedure for resolv- 
ing the experimental lifetime spectrum into three or 
me re components is described. Previous attempts 
to develop procedures and programs for analyzing 
mi Iticomponent positron lifetime spectra are de- 
scribed in references 32 to 34. 

Sp ectrum Analysis 

Procedure 

The analysis procedure is described by refer- 
ring to the experimental data shown in figure 3(b). 
Thu data have been plotted on a semilog paper in 
orcer to simplify the discussion. A cursory exam- 
ination of this figure shows that it has three dis- 
tin t components superimposed on a constant back- 
ground. The stripping process starts by subtracting 
the background, which can be easily determined by 
av< raging about 50 channels to the far right where 
no genuine positron decay events contribute to the 
observed data. Figure 4(a) shows the spectrum re- 
maining after the background is subtracted. (The 
tin e-scale origin t z has been shifted to channel zero 
and the abscissa has been expanded for the sake of 
cla ity and for convenience of discussion.) Attention 
is irst focused on the longest lifetime component 
(i.e , the third component). The first and second 
components do not seem to make any contributions 
to he data beyond channel 150. It would therefore 
appear that data in channels 150 to 180 are solely 
attributable to the third component, the initial es- 
timate of which can be obtained by a least-squares 
fit to the data in these channels. After subtract- 
ing the third component from the background-free 
spectrum, the remaining spectrum (fig. 4(b)) has 
twe distinct components. The counts in channels 
40 to 60 are reasonably free from the effects of the 
firs component. A least-squares fit to the data in 
channels 40 to 60 then provides the necessary ini- 
tial information about the second component. Fig- 
ure 4(c) illustrates the fit for the second component. 
Th spectrum remaining after subtracting the sec- 
onc component is then solely attributable to the first 
component. Figure 4(d) shows this residual spec- 
tru n. Analysis of the residual spectrum then pro- 
vides an initial estimate for the first component in the 
spe drum. Figure 4(e) illustrates the fit for the first 
component. After obtaining initial estimates for the 
thme components, the finite- resolution effects are in- 
cluded in each component as discussed in the preced- 
ing section. The sum of these modified components 
gives the initial computed least-squares spectrum 
for comparison with the experimental spectrum. 


3 



Figure 4(f) illustrates the initial computed spectrum. 
This initial computed spectrum is then reiterated un- 
til the difference between the computed and the ex- 
perimental spectra is minimized. The final compar- 
ison between the fitted and experimental spectra is 
shown in figure 4(g). A computer program devel- 
oped by using these procedures is described in ap- 
pendixes A and B. 

Test of New Method 

To validate the computational procedure de- 
scribed in appendixes A and B, several artificial 
sets of data with random fluctuations were first con- 
structed. An experimentally observed Co 60 spectrum 
was used to convolve the input sum of exponentials 
for simulation of the true experimental data. These 
simulated spectra were then analyzed using the same 
procedure as for the actual data. A second experi- 
mental Co 60 spectrum was used as the input for the 
resolution function for PAPLS and POSITRONFIT 
analyses. Typical results for two such cases are sum- 
marized in tables I to III. 

It is apparent from the results summarized in ta- 
bles II and III that the present method deconvolves 
the artificial spectra correctly. It is also apparent 
that PAPLS provides closer fits to the input values 
for almost all the parameters than the POSITRON- 
FIT program, particularly in the case of intensity val- 
ues for the three components. Figures 5(a) and 5(b) 
illustrate the comparison between the fits predicted 
by PAPLS and POSITRONFIT programs, respec- 
tively, for case 1. 

Applications 

The PAPLS program has been applied in the 
analysis of positron lifetime spectra observed in 
metal-ion containing Epon 828 epoxy samples. The 
target specifications and experimental lifetime spec- 
tra observed in them are listed in appendix C. 
The same spectra were also analyzed using the 
POSITRONFIT program for the purpose of com- 
parison. The results of this analysis are also in- 
cluded in table IV. A comparison of the two sets 
of values shows that while they are in general 
agreement, certain differences do exist. For exam- 
ple, most of the POSITRONFIT lifetime values are 
slightly higher than the values given by the PAPLS 
program. 1 Similarly, the intensities of all the 


1 The PAPLS program gives slightly higher r<i and 
values in Epon 828 containing 0.1 mole fraction of 
Cr (DMSO) 6 (C10 4 ) 3 . 


POSITRONFIT shortest lifetime components are 
higher, whereas the intensities of the intermediate 
lifetime components are lower than the correspond- 
ing PAPLS values. These differences are attributable 
to the fact that the PAPLS program uses the ex- 
act resolution function, whereas the POSITRONFIT 
program assumes it to be a Gaussian or sum of 
Gaussians. The resolution function is not rigorously 
representable by a sum of Gaussians or Gaussians 
and exponentials. 

Because of the manner in which t z is calcu- 
lated in PAPLS, it takes considerably longer to 
complete a PAPLS analysis than a POSITRONFIT 
analysis. Overall, however, the PAPLS program 
is more accurate, because it uses an exact resolu- 
tion function rather than an approximate analytical 
representation. 

Concluding Remarks 

A new technique for analyzing multicomponent 
positron lifetime spectra in polymers has been de- 
veloped. It utilizes the actual experimental timing- 
system- resolution function, rather than an analyti- 
cal approximation, for deconvolving the experimental 
spectra. These concepts have been incorporated into 
a computer program for analyzing positron lifetime 
spectra (PAPLS) in polymers and other condensed 
media where positronium is likely to be formed. The 
validity of this program has been tested using sev- 
eral artificially generated data sets with random sta- 
tistical fluctuations. Typically, the fitted parame- 
ters agree with the input values within ±2 percent. 
These same data sets were also analyzed using the 
widely used POSITRONFIT program for the pur- 
pose of comparison. In almost all cases, the present 
technique gives closer fits to the input values. Both 
PAPLS and POSITRONFIT have been used to an- 
alyze several lifetime spectra measured in metal-ion 
containing Epon 828 epoxy samples. Even though 
the results obtained by the two programs are essen- 
tially equal, the PAPLS results are considered more 
accurate, because PAPLS does not make an a pri- 
ori assumption about the nature of the life timing- 
system-resolution function. 

NASA Langley Research Center 
Hampton, Virginia 23665-5225 
October 11, 1988 
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Appendix A 

Computer Program (PAPLS) 

The computer program for analyzing positron 
lifetime spectra (PAPLS) is written in FORTRAN 
Version 5 language for the Control Data CYBER 170 
series digital computer system with the Network Op- 
erating System (NOS 2.4). Machine dependence is 
limited to the use of several library routines (matrix 
inversion and graphics routines), and the program 
should be readily adaptable to other computer sys- 
tems. All routines not included in the program are 
described in the code. The program requires approx- 
imately 70000 octal locations of core storage, and a 
typical analysis requires approximately 100 central 
processing unit (CPU) seconds. 

It is assumed in the program that the positron 
lifetime spectrum can be described by the equation 


n k 


L?=0 t=l 


+ B (Al) 


where n is the number of counts in channel Ar, 
is the amplitude at time zero t z of each component, 
Aj = l/r^, t x is the lifetime of each component, tj is 
the time corresponding to channel j, B is a constant 
background, and the values of W j ^ are normalized 
weights representing the resolution function smear- 
ing between channels j and k. In this expression, the 
values of the amplitudes A x are zero for t < t z . The 
number of components is arbitrary, but the particu- 
lar version 'of the program being described has three 
components. The weights may be defined by us- 
ing either an experimental resolution function or an 
analytical representation. The program defaults to a 
Gaussian time function if an experimental spectrum 
is not provided. The default weights are computed 
by 

W ik = -^=e-^- t ^ 2 /° 2 (A2) 

(Ty/TT 

where a is the standard deviation. 

The process of obtaining a solution begins by first 
eliminating the background from the spectrum. The 
background is obtained by averaging counts over sev- 
eral channels (the default is 51 channels) at the high 
end of the spectrum. Thus, it is important that 
enough channels are provided so that the contribu- 
tion of genuine positron decay events is negligible 
within the region used for background averaging. 

Before applying the least-squares technique, de- 
scribed in detail in appendix B, to equation (Al) with 
B subtracted from both sides of the expression, it is 
necessary to estimate initial values for A x and X t . 
Rather crude estimates are usually sufficient when 


he experimental data are relatively noise free. A 
stripping process is used to obtain these estimates. 
For this purpose, equation (Al) is approximated by 
he following simplistic expression: 

In (n - B) = ln(Ai) - A it + ln(A 2 ) - A 2 / 

+ In(A 3 )-A 3 £ (A3) 

\t this point of the solution process, it is further 
issumed that the resolution is perfect and that the 
^ is exactly at the center of a channel a fixed number 
>f channels (the default is 2 channels) to the left of 
he observed spectrum maximum. 

The longest lifetime component is first addressed 
>y examining only the data in a range where the 
contribution from the other two components has 
decayed to a negligible value (the default range is 
rom channel 100 to channel 130 to the right of t z ). 
The data in these channels are least-squares fitted 
o a straight line to provide the initial estimates for 
I 3 and A 3 . Using these values, the third component 
3 subtracted from equation (A3). The intermediate 
ifetime component is next considered by examining 
• »nly the data in a range where the short lifetime 
- omponent has decayed to a small value (the default 
ange is from channel 38 to channel 58 to the right 
of t z ). As before, a straight line fit yields the initial 
estimates for A 2 and A 2 . The intermediate lifetime 
component is next subtracted from equation (A3), 
finally, initial estimates for A\ and Ai are obtained 
by a straight line fit to the remaining data. The 
default range for this fit is from channel 11 through 

< hannel 20 to the right of t z . 

With these six values used as initial estimates, an 
iterative process is initiated to determine the solution 
of equation (Al). The first step is to hold these six 
' alues constant and find the value of t 2 , which min- 
imizes the sum of the squares of the residuals. The 

< ntire range of data is used starting eight channels 
to the left of the peak. The initial estimates for A x 
i nd A^ are then used with this value of t z , and the 
bast-squares technique is employed to generate new 
values for A x and A t . For this portion of the solu- 
tion, the range of data used does not include data 
irom eight channels to the left of the peak to three 
( hannels to the right of the peak. The improved val- 
ves of A x and \ x are then held fixed, and t z is once 
1 lore adjusted to minimize the sum of the squares of 
1 he residuals. This iterative process continues until 
i\t z < 10 ^ channel, where A t z is the change in t z 
Irom one iteration to the next. 

The resulting seven values of Aj, and t z are 
t hen used as initial estimates for one final iterative 
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procedure. This procedure is identical to the pre- 
vious step, except that all the data are used when 
employing the least-squares technique. 

Program Input 

All program control input is provided through 
the NAMELIST name INDAT. A summary of the 
input variables is provided below. (Default values 
are provided in parentheses.) 

NDIV (51): 

Number of channels at high end of spectrum to 
be used in obtaining average background. 

BKG (computed on basis of NDIV): 

Average background. If specified, it will override 
the computed value. 

TSF (61.27): 

Time-scale factor in psec per channel. 

NPO (120,38,11): 

Array that provides channel offsets for defining 
regions used in obtaining initial estimates. 

NPL (30,20,9): 

Array that provides interval lengths for defining 
regions used in obtaining initial estimates. De- 
faults are based on the default value of TSF. 

IDCHAN (-2): 

Correction factor to compensate for apparent 
shift in peak location that results from time- 
resolution function. In general, the true value 
of t z is to the left of the peak in the experimen- 
tal data. An examination of the program output 
should verify that the computed and experimen- 
tal data peak in the same channel; if not, an ap- 
propriate value of IDCHAN will correct for this. 

IGAUSS (1): 

Flag to indicate if a Gaussian time-resolution 
function is to be used in fit. If IGAUSS = 0, 
no time resolution function is used. 

SIGINV (0.135621): 

Value of 1 /ct 2 , where a is standard deviation of 
Gaussian time-resolution function. 

ITERMX (30): 

Maximum number of iterations to be allowed in 
least-squares fit. 

IPRFLG (0,0, 1,0, 0,1): 

Array controlling quantity of printed output. 
Each of six flags controls a specific output item. 
If the flag is “1,” the item is printed. If the flag is 
“0,” the item is not printed. The flags correspond 
to the following output items: 

1. Raw data 

2. Background corrected data 


3. Initial-estimate fit results 

4. Spectrum with initial fit subtracted 

5. Fit using initial estimates 

6. Final fit 

IPLFLG (1,0, 1,1, 1,1): 

Array that controls the quantity of graphic out- 
put. Each of six flags controls a specific plot. If 
the flag is “1,” the plot is generated. If the flag 
is “0,” the plot is not generated. The flags corre- 
spond to the same items as those for IPRFLG. 

IPEEK (0,0): 

Array used to simply examine data. If 
IPEEK(l) / 0, the program produces a plot of 
the log of the data versus channel number; the 
plot begins with channel number IPEEK(l) to the 
right of the peak and continues through channel 
number IPEEK(2) to the right of the peak. No fit 
is attempted. This feature is useful for determin- 
ing the ranges used for obtaining initial estimates. 

In addition to program control input, the pro- 
gram reads the spectrum from TAPE1. The first 
record on this file should contain the number of chan- 
nels and a data identification label in 14, A7 format. 
No more than 800 channels are allowed. The re- 
maining records on this file contain integer counts 
for the number of channels specified in list-directed 
format. When an experimental resolution function 
spectrum is provided, these data complete the con- 
tents of TAPEl, which is also in list-directed format. 

Program Output 

The bulk of the program output is controlled by 
user-defined flags as indicated in the section “Pro- 
gram Input.” The initial output (fig. Al) is always 
produced and provides an echo of the input spec- 
trum with the data identification label, the number 
of channels, the value of the average background, and 
the number of channels used to compute the average 
background. If desired, the raw data are then out- 
put, as illustrated in figure A2, and give the time 
in psec, the counts, and the log of the counts for 
each channel. The data that are corrected for back- 
ground may also be printed in the same four-column 
format. Data relating to the spectrum peak are then 
output, as shown in figure A3, and give the max- 
imum value of the background-corrected spectrum, 
the channel where the peak occurs, and the channel 
corresponding to t z . For each component, a sum- 
mary of the initial-estimate calculations can be pro- 
duced as shown in figure A4. This summary identifies 
the component, the channel range used to obtain the 
estimates, and the values (with errors in parenthe- 
ses) for the slope A, lifetime TAU, intercept ln(A), 
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and the area under the resulting curve. Also output 
are the fit standard deviation, the number of itera- 
tions required for convergence, and an error code. If 
this error code is not zero, the solution did not con- 
verge. For each component, the spectrum that re- 
sults from subtracting the initial-estimate curve can 
also be output in a form similar to that in figure A2. 
After initial estimates have been determined for all 
six unknowns, the spectrum calculated from these 
values can be output for comparison with the exper- 
imental data. A sample of this output is included 
as figure A5; the sample provides the channel num- 
ber, time (in psec), experimental data, and the fit- 
ted spectrum by using initial estimates. These data 
can be examined to verify the experimental data and 
the fit peak in the same channel. If not, the input 
variable IDCHAN requires adjustment. Sample out- 
put corresponding to the final-fit values is illustrated 


in figure A6. For each component, values are given 
(with errors in parentheses) for the slope A, life- 
time TAU, intensity ^Area; j J2 Area^ , intercept 

ln(A), and the area under the resulting curve. Values 
are also provided for the fit standard deviation, the 
number of iterations required for convergence, and an 
error code. The final output is similar to figure A5, 
where the calculated spectrum is obtained using the 
final- fit values. 

In addition to the printed output, the program 
produces graphic output, all of which is controlled 
by user-defined flags. The plots corresponding to 
the six flags are shown in figure A7. Results similar 
to figures A7(c) and A7(d) are produced for each 
component. 
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Program Listing 


PROGRAM PAPLS(OUTPUT , INPUT ,TAPE5= IN PUT ,TAPE6=OUTPUT , 

I TAPE1) 

C******* ******************************************************* ********* 


C * 
C PAPLS IS A COMPUTER PROGRAM FOR ANALYZING POSITRON * 
C LIFETIME SPECTRA. THE CURRENT VERSION OF THE PROGRAM * 
C IS LIMITED TO THREE COMPONENTS, BUT MINOR CODING * 
C CHANGES WOULD BE REQUIRED TO ENABLE THE PROGRAM * 
C TO PROCESS MORE OR FEWER COMPONENTS. AN ADVANTAGE * 
C OF PAPLS OVER OTHER PROGRAMS OF A SIMILAR NATURE * 
C IS ITS HANDLING OF THE TIME RESOLUTION FUNCTION. * 
C THE RESOLUTION FUNCTION IS DESCRIBED IN TABULAR * 
C FORM, THUS ALLOWING ANY ANALYTIC (OR NON-ANALYTIC) * 
C REPRESENTATION. FOR A NON-ANALYTIC FUNCTION, A LOW * 
C LEVEL OF NOISE IN THE REPRESENTATION OF THE TIME * 
C RESOLUTION IS REQUIRED FOR RELIABLE PROGRAM * 
C EXECUTION. * 
C * 


C *********************************************************************** 
COMMON/INDAT/ICNT(800) ,NPTS , NPTSS , BKG , NDI V , IMAXC ,TSF , IDCHAN , 

I IPRFLG(LO) ,IPLFLG(IO) ,NPO(4) ,NPL(4) , LABEL, CFRAC 
COMMON/PEEK/I PEEK( 2) 

COMMON/GAUSDAT/AMP , SIGINV , IGAUSS ,GAUSS(9 I ) 

COMMON /OTD AT /SLOPEX(4 ) ,XOFFX(4) , ESLOPE(4) , EXOFF(4) , 

I NP1ST(4) ,NPLST(4) ,TAU(4) , ETAU(4) ,AREA(4) ,EAREA(4) 
COMMON/F'ITDAT/XX(800) ,YY(800) 

COMMON/PRDAT/DATAX(800) ,DATAY(800) ,DATAXT(800) ,DATAYL(800) 
COMMON/CLCDAT/X(800) ,Y(800) ,YL(800) 

COMMON/CHARS /MSGPR( 4 ) 

DIMENSION SOL(7) ,E(7) 

COMMON/NEWCOM/SOL,E ,ITER2 ,IERR,STD 
DIMENSION Z(800) 

DATA MSGPR/6H FIRST , 6HSECOND , 6H THIRD , 6HFOURTH/ 

DATA SLOPEX ,XOFF’X/ 8*0 . / 

DATA NPASS/O/ 

1 FORMAT ( 1H1 , 5X, ' INPUT SPECTRUM (',A7 ,')'//) 

2 F0RMAT(5I 10) 

3 FORMAT (LH0,5X, 'NUMBER OF CHANNELS IN DATA FILE = ',15) 

4 FORMAT (1 HO, 5 X,' BACKGROUND AVG. = ' , IPEI3 . 5 , 5X, 'NDIV = ',13) 

5 F0RMAT(1H0,5X, 'NUMBER OF COMPONENTS TO FIT = ',15) 

6 FORMAT ( 1HI , 5X, 'RAW SPECTRUM'// 

I 5X, 'CHANNEL' , IOX, 'TIME(PS) ' , 10X, 'COUNTS' , 10X, 'LOG(COUNTS) ' // ) 

7 FORMAT (IX ,F10.0,7X,LPEi3.5,5X,lPE13.5,5X,lPEL3.5) 

8 F0RMAT(1H1,5X, 'BACKGROUND CORRECTED SPECTRUM'// 

1 5X, 'CHANNEL' , IOX, 'TIME(PS) ' , IOX, 'COUNTS' , IOX, 'LOG(COUNTS) '// ) 

9 F0RMAT(1HI,5X, 'SPECTRUM MAXIMUM = ',1PE13.5/ 

1 1 HO, 5X, 'CHANNEL NUMBER CORRESPONDING TO SPECTRUM MAXIMUM = ', 

2 15/ 1H0.5X, 'CORRECTED CHANNEL = ',15) 

10 FORMAT (1 HO, 5X, 'SPECTRUM MAXIMUM = '.1PE13.5/ 

1 1H0,5X, 'CHANNEL NUMBER CORRESPONDING TO SPECTRUM MAXIMUM = ', 

2 I5/1H0,5X, 'CORRECTED CHANNEL = ',15) 

11 F0RMAT(1H1,5X, 'SPECTRUM CORRECTED FOR ',A6,' FIT'// 

1 5X, 'CHANNEL' , IOX, 'TIME (PS)' , IOX, 'COUNTS' , IOX, 'LOG(COUNTS)'// ) 
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12 FORMAT ( mi ,5X,A6 , ' FIT (FROM', 15,' Ti ' , 15 , ' ) '/// 

1 10X, 'SLOPE = ' , 1PE16 . 8 , IX, 

2 '(',1PE16.8,')'//10X,'TAU = ' , 1PE16 .8, IX, ' (' , 

3 1PE16.8,')'//10X, 'INTERCEPT = ' , 1PE16 . 8 , IX, ' ( ' , 

4 1PEL6.8,')'//10X, 'AREA = ' , 1PE16.8 , IX, ' ( ' , 

5 1PEL6.8,')'//10X, 'STANDARD DEVIATION = '.1PE16.8// 

6 10X, 'ITERATIONS = ' .I5//10X, 'IERR = '15///) 

13 FORMAT (1 til, 5X,' FITTED SPECTRUM'// 

1 5X, 'CHANNEL' ,1 OX, 'TIME (PS)', 1 OX, 'COUNT S' , 10X, ' FIT '//) 

14 F0RMAT(1H1,5X, 'NONLINEAR FIT'/// 

1 10X,' SLOPE = ' ,3( 1PE14. 6 , IX, ' (+' , 1 PE14 .6 , ') ' ,3X) / 

A 1H+,21X,3(16X,'_' , 18X)/ 

2 10X,' TAU = ',3(1 PE 14. 6, IX, '(+',] PE 14. 6,')' ,3X) / 

A 1H+,21X,3(16X, '_' , 18X)/ 

3 10X, 'INTERCEPT = ' , 3( 1PE14 .6 , IX, ' (+' , 1 PE14. 6 , ' ) ' ,3X) / 

A 1H+, 21X, 3( 16X, '_' , 18X)/ 

4 10X,' AREA = ', 3(1 PE 14. 6, IX, '(+',1 PE 14. 6,')', 3X)/ 

A 1H+,21X,3(16X,'_',18X)/ 

5 10X, 'INTENSITY = ' , 3( 1PE14 . 6 , IX , ' (+' , 1 PE14 . 6 , ' ) ' , 3X) / 

A 1H+,21X,3(16X, ' ' ,18X)/ 

6 10X,' TO = ', 1 PE 14. 6, IX, '(+',1 PE 14. 6,')'/ 

A 1H+,21X,1(16X,'_',18X)/ 

7 10X, 'STANDARD DEVIAIION = '.1PE14.6// 

8 10X, 'ITERATIONS = ' , 15 , ' , ' , 15/ / 

8 10X,' ICT = ',15// 

9 10X, 'IERR = ',15///) 

15 FORMAT(1HO,5X, 'LINEAR FIT STANDARD DEVIATION = '.1PE16.8, 

1 5X, 'VARIANCE = '.1PE16.8) 

16 FORMAT (1H0.5X, 'NONLINEAR FIT STANDARD IEVIATION = ',1PE16.8, 

1 5X, 'VARIANCE = '.1PE16.8) 

17 FORMAT (1H0,5X, 'TOTAL COUNTS = ',110) 

18 FORMAT ( ' GTOT = ',Elb.8) 

c 

C 

C SUBROUTINE PSEUDO INITIALIZES THE 1 LOT VECTOR FILE. 

C THIS SUBROUTINE IS DOCUMENTED IN "LANGLEY GRAPHICS 

C SYSTEM", CENTRAL SCIENTIFIC COMPUTING COMPLEX 

C DOCUMENT G-3 

C 

C 

CALL PSEUDO 

C-H-H-+ 

C 

READ INPUT DATA 

C l- EL EE 

ICT = 0 

CALL DATRD 

WRITE (6,1) LABEL 

WRITE (6, 2) (ICNT(I) ,1=1 ,NPTS) 

ISUMCT = 0. 

DO 20 1=1 , NPTS 

ISUMCT = ISUMCT + FLOAT(ICNT(I) ) 

20 CONTINUE 

WRITE(6 , 17 ) ISUMCT 
WRITE( 6 ,3) NPTS 
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WRITE ( 6 ,4) BKG.NDIV 
NFIT = 0 

DO 30 1=1,4 

IF(NPL(I).NE.O) NFIT = NFIT + 1 
30 CONTINUE 

WRITE(6 , 5) NFIT 

C-H-H-+ 

C 

C STORE SPECTRUM IN DATAY 

C STORE LOG( SPECTRUM) IN DATAYL 

C STORE CHANNEL NUMBER IN DATAX 

C STORE TIME (PS) IN DAT AXT 

C 

CH-H4- 

CNTTOT = 0. 

DO 40 1=1 , NPTS 

CNTTOT = CNTTOT + ICNT(I) 

DATAY ( I ) = FLOAT ( I CNT ( I ) ) 

DATAYL(I) = 0. 

IF ( DATAY ( I) .GE . 1 . ) DATAYL ( I ) = ALOG(DATAY( I ) ) 
DATAX(I) = FLOAT ( I ) 

DAT AXT ( I ) = DATAX( I ) *TSF 
40 CONTINUE 

c um- 

c 

C WRITE OUT RAW DATA IF REQUESTED 

C 

Cl 1 PH - 

IF ( IPRFLG( 1 ) .EQ . 0) GO TO 60 
WRITE(6 , 6) 

DO 50 1=1 ,NPTS 

WRITE ( 6 , 7 ) DATAX ( I ) , DATAXT ( I) , DATAY (I) , DATAYL ( I ) 
50 CONTINUE 
60 CONTINUE 
C HID 
C 

C PLOT RAW SPECTRUM IF REQUESTED 

C 

G H l i t - 

IF ( IPLFLG( 1 ) . EQ. 0) GO TO 70 

CALL DATPLT (DATAX , DATAYL , 1 , SLOPE , XOFF , NPL ( 1 ) , 

1 NPO(l).MFIT) 

70 CONTINUE 

IF(IPEEK( 1) . NE .0) GO TO 400 

C I I- l- H - 

C 

C ELIMINATE BACKGROUND 

C 

CH+H 

DO 80 1=1 ,NPTS 

DATAY ( I ) = DATAY ( I ) - BKG 

DATAYL(I) = 0. 

IF(DATAY( I ) .GE . 1 . ) DATAYL(I) = ALOG(DATAY( I ) ) 

80 CONTINUE 
H -t- H - 
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C WRITE OUT BACKGROUND CORRECTED DAT/ IF REQUESTED 

C 

C-H-H-+ 

IF(IPRFLG(2) .EQ.O) GO TO 100 
WRITE( 6 ,8) 

DO 90 1=1 , NPTS 

WRITE ( 6 , 7 ) DATAX( I ) , DATAXT ( I ) ,DATAY( I ) , DATAYL( I ) 

90 CONTINUE 
100 CONTINUE 
■H-+-H- 


PLOT BACKGROUND CORRECTED DATA IF REQUESTED 
C 

C++-R++ 

IF(IPLFLG(2) .EQ.O) GO TO 110 

CALL DATPLT (DAT AX , DATAYL , 2 , SLOPE , XOFF , N PL ( 1 ) , 

1 NPO( 1 ) , MFIT ) 

110 CONTINUE 
CR H I I 
C 

C LOCATE CHANNEL CONTAINING SPECTRUM MAXIMUM 

C 

C-H-R++ 

YMAX = 0. 

DO 120 1=1 ,NPTS 

IF(DATAY(I) .LT.YMAX) GO TO 120 
YMAX = DATAY(I) 

IMAX = I 
120 CONTINUE 
C R R R RR 
C 

C SET CORRECTED CHANNEL AND SHIFT DAT \ 

C STORE DATA IN CALCULATION ARRAYS 

C 

C R R + -R + 

IMAXC = IMAX + I DC HAN 
NPTSS = NPTS - IMAXC + 1 
DO 130 1 = 1, NPTSS 
X(I) = FLOAT ( I- 1) 

Y(I) = DATAY ( I+IMAXC- 1 ) 

YL( I ) = DATAYL (I+IMAXC-1) 

130 CONTINUE 

IF(IPRFLG(1)+IPRFLG(2) .EQ.O) WRITE(6,10) YMAX , IMAX , IMAXC 
IF(IPRFLG(l)+IPRFLG(2).NE.O) WRITE (6, 9) YMAX , IMAX , IMAXC 

C RRRRR 

C 

C CALCULATE GAUSS WEIGHTS 

C 

CRRRRR 

ICOR = IMAXC - 46 
GTOT = 0. 

IF(IGAUSS.EQ.2) GO TO 150 
DO 140 1=38,54 
XI = FLOAT ( I+ICOR) 

GAUSS (I) = AMP*EXP(-SIGINV*(XI - FLOAT( LMAXC))**2) 



nooo ooooo nooo 


GTOT = GTOT + GAUSS(I) 

140 CONTINUE 
GO TO 170 
150 CONTINUE 

DO 160 1=38,54 

GTOT = GTOT + GAUSS (I) 

160 CONTINUE 
GO TO 190 
170 CONTINUE 

IF( IGAUSS . EQ . 1 ) GO TO 190 
DO 180 1=1,91 
GAUSS(I) = 0. 

180 CONTINUE 

GAUSS(46) = 1. 

190 CONTINUE 

WRITE(6 , 18) GTOT 

C-H--H-+ 


LOOP THROUGH DESIRED FITS 


T-H-H- 

MFIT = 0 
DO 240 K=1 ,4 

IF(NPL(K.) .EQ.O) GO TO 240 
MFIT = MFIT + 1 

CALL LSQ2(NPL(K) ,X(NP0(K) ) ,YL(NP0(K.) ), SLOPE ,XOFF , 
1 ITER, IERR,STD ,E1 , E2) 

SLOPEX(MFIT) = SLOPE 
XOFFX(MFIT) = XOFF 
ESLOPE(MFIT) = El 
EXOFF(MFIT) = E2 
TAU(MFIT) = -1. /SLOPE 
TAU(MFIT) = TAU(MFIT)*TSF 
ETAU(MFIT) = El/SLOPE**2 
ETAU(MFIT) = ETAU(MFIT)*TSF 
AREA(MFIT) = TAU(MFIT)*EXP(XOFF) 

EAREA(MFIT) = ETAU(MFIT)*EXP(XOFF) 

1 + EXOFF(MFIT)*TAU(MFIT) 

NPIST(MFIT) = NPO(K) 

NPLST(MFIT) = NPO(R) + NPL(K) 

I H H - 


WRITE OUT FIT RESULTS IF REQUESTED 


HH -+ 

IF(IPRFLG(3) .EQ.O) GO TO 200 

WRITE(6 , 12) MSGPR(MFIT) , NPIST(MFIT) , NPLST(MFIT), SLOPE, El, 

1 TAU(MFIT) .ETAU(MFIT) ,XOFF,E2 , AREA(MFIT) ,EAREA(MFIT) ,STD, 

2 ITER,IERR 
200 CONTINUE 

H - H -t - 


PLOT THE FIT IF REQUESTED 


C+++++ 
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IF( IPLFLG( 3 ) . NE .0) CALL DATPLT(X , YL, 5 SLOPE, 

1 XOFF ,NPL(K) ,NPO(K) ,MFIT) 

C+-H-H- 

C 

C SUBTRACT OUT THIS RESULT 

C 

Cl H H 

FN = EXP(XOFF) 

DO 2 LO 1=1 , NPTSS 

Y( I ) = Y( I ) - FN*EXP(SLOPE*X( I ) ) 

I F ( Y ( I ) . LT . 0 . ) Y( I ) = 0. 

YL( I ) = 0. 

IF ( Y( I ) .GT . 1 . ) YL( I ) = ALOG( Y( I ) ) 

210 CONTINUE 
C-H-+++ 

C 

WRITE OUT CORRECTED SPECTRUM IF REQUESTED 


l -H - H 

IF( IPRFLG(4 ) . EQ • 0) GO TO 230 
WRITE (6,11) MSGPR(K) 

DO 220 1=1, NPTSS 
IF(I.LT.IMAXC) 

1WRITE(6, 7) DATAX(I) ,DATAXT(I) ,DATAY(i; ,DATAYL(I ) 
J = I - IMAXC + 1 
IF(I.GE.IMAXC) 

1WRITE(6 , 7 ) DATAX(I) ,DATAXT(I) ,Y( J) ,YL( J) 

220 CONTINUE 

-H-H-+ 


PLOT THE CORRECTED SPECTRUM IF REOUESTED 


1 H - H 

IF(IPLFLG(4) .EQ.O) GO TO 230 

CALL DATPLT(X,YL, 3, SLOPE, XOFF, NPL(K) , 

1 NPO(K) ,MFIT) 

230 CONTINUE 
240 CONTINUE 
250 CONTINUE 

NPASS = NPASS + 1 

I I I I I- 

RESTORE DATA 

H-H-+ 

DO 260 1=1 , NPTS 
DATAY(I) = DATAY(I) + BKG 
DATAYL ( I ) = 0. 

IF(DATAY(I).GE.l.) DATAYL(I) = AL0G(D7 TAY( I ) ) 
260 CONTINUE 

l -H - H- 


COMPUTE FITTED SPECTRUM 
C 

C+ l- H I 
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DO 290 1=1 ,NPTS 
ZZZ = 0. 

YY( I) = 0. 

XX(I) = FLOAT ( I ) 

IF ( I . LT . IMAXC ) GO TO 280 

ARGX = DATAX(I) - ( FLOAT ( IMAXC) + CFRAC ) 

IF(I.EQ. IMAXC) ARGX = 0.5*(0.5 - CFRAC) 

I F ( ARGX . LT . 0 . ) ARGX = 0. 

DO 270 KK=1 ,MF IT 

ZZZ = ZZZ + EXP(XOFFX(KK))*EXP(SLOPEX(KK)*ARGX) 

270 CONTINUE 

IF( I . EQ. IMAXC ) ZZZ = (0.5 - CFRAC)*ZZZ 
280 CONTINUE 

ZZZ = ZZZ + BKG 

IF (ZZZ . LE . 1 . ) GO TO 290 

YY( I ) = ALOG(ZZZ) 

290 CONTINUE 

IF(IGAUSS.EQ.O) GO TO 330 
IBEG = IMAXC - 8 
IEND = NPTS-10 
DO 310 I=IBEG, IEND 
GCOR = 0. 

DO 300 K=38,54 
INDX = 1+4 6 -K 
IF(INDX.LE.O) GO TO 300 
YDATA = EXP (YY( INDX)) - BKG 
GCOR = GCOR + GAUSS (K)*YDATA 
300 CONTINUE 

Z(I) = GCOR + BKG 
Z(I) = ALOG ( Z ( I ) ) 

310 CONTINUE 

DO 320 I=IBEG, IEND 
YY( I ) = Z ( I ) 

320 CONTINUE 
330 CONTINUE 

STDTMP = 0. 

VARTMP = 0. 

NBEG = IMAXC - 6 
DO 340 I=NBEG, NPTS 

STDTMP = STDTMP + (EXP(YY(I)) - FLOAT( ICNT( I) ) )**2 
OBSRV = 1. 

IF( ICNT ( I ) . NE .0) OBSRV = 1 . /FLOAT( ICNT( I ) ) 

VARTMP = VARTMP + OBSRV*( EXP( YY ( I) ) - FLOAT (ICNT(I)))**2 
340 CONTINUE 

STDTMP = SQRT( STDTMP /FLOAT (NPTS-NBEG) ) 

VARTMP = VARTMP/ (FL0AT(NPTS-NBEG+1 - 7)) 

I F(NPASS . EQ . 1 ) STDLIN = STDTMP 
IF (NPASS . EQ . 1 ) VARLIN = VARTMP 
IF(NPASS . EQ . 2 ) STDNLIN = STDTMP 
IF(NPASS . EQ . 2) VARNLIN = VARTMP 

C+++++ 

C 

C WRITE OUT THE FITTED SPECTRUM IF REQUESTED 

C 

C+++++ 
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IF(IPRFLG(NPASS+4) .EQ.O) GO TO 360 
WRITE(6 ,13) 

DO 350 1=1 ,NPTS 
FIT = EXP(YY(I) ) 

WRITE(6 , 7 ) DATAX( I ) , DATAXT ( I ) , DATAY ( I ' ,FIT 
350 CONTINUE 
360 CONTINUE 

FITMAX = EXP(YY(IMAXC) ) 

RELMAX = DATAY (IMAXC) 

RATMAX = RELMAX/ FITMAX 

IF(NPASS .EQ. 1) WRITE(6 , 15) STDLIN , VARI IN 
IF(NPASS . EQ. 2) WRITE(6,16) STDNLIN.VAKNLIN 

C l i nt - 

C 

PLOT THE FITTED SPECTRUM IF REQUESTED 


H-H -t- 

IF(IPLFLG(NPASS+4) .EQ.O) GO TO 370 

CALL DATPLT(DATAX,DATAYL,4*NPASS , SLOPS ,XOFF , NPL( 1 ) , 
1 NPO(l).MFIT) 

370 CONTINUE 

IF(NPASS . EQ . 2) GO TO 400 

I I I H- 

DETERMINE NONLINEAR FIT 

C+-H-H- 

SOL(l) = EXP(XOFFX( 1) ) 

S0L(1) = S0L( 1)*RATMAX 
SOL(2) = -SLOPEX(l) 

SOL(3) = EXP(XOFFX(2) ) 

SOL(3) = S0L(3)*RATMAX 
SOL(4) = -SLOPEX(2) 

SOL(5) = EXP(XOFFX(3) ) 

S0L(5) = SOL( 5)*RATMAX 
SOL(6) = -SL0PEX(3) 

SOL(7 ) = 0.5 
DO 380 I=L,NPTS 
DATAY ( I ) = DATAY ( I ) - BK.G 
380 CONTINUE 

SOLI = SOL(l) 

SOL2 = S0L(2) 

S0L3 = SOL(3) 

SOL4 = SOL(4) 

SOL5 = SOL(5) 

SOL6 = SOL (6) 

SOL7 = S0L(7) 

390 CONTINUE 

SOL(l) = SOLI 
SOL(2) = S0L2 
SOL(3) = S0L3 
S0L(4) = S0L4 
SOL(5) = S0L5 
SOL(6) = SOL6 
S0L( 7 ) = S0L7 
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CALL LSQ6 ( NPTS , X , DATAY , SOL , E , ITER1 , 1 ERR , STD , IMAXC , 
1 NPO,NPL, CFRAC ,ICT) 

SOLL = S0L(1) 

SOL2 = SOL( 2) 

SOL3 = SOL( 3) 

SOLA = SOL(4) 

SOL5 = SOL( 5) 

SOL6 = S0L(6) 

S0L7 = SOL(7 ) 

SLOPEX(l) = - SOL(2) 

SLOPEX(2) = - SOL(4) 

SLOPEX( 3) = - SOL(6) 

ESL1 = E( 2) 

ESL2 = E(4) 

ESL3 = E(6) 

EXOl = ALOG( E ( 1 ) ) 

EX02 = ALOG(E(3) ) 

EX03 = ALOG(E(5) ) 

SOL(2) = S0L(3) 

SOL(3) = SOL(5) 

SOL(4) = CFRAC 
SOL(5) = -SLOPEX(l) 

SOL(6) = -SLOPEX(2) 

SOL(7) = -SLOPEX( 3) 

CFRACO = CFRAC 
CALL LSQ1 

XOFFX(l) = ALOG( SOL( 1 ) ) 

XOFFX(2) = ALOG(SOL(2) ) 

XOFFX(3) = ALOG(SOL(3) ) 

TZERO = FLOAT (IMAXC) + SOL(4) +0.5 
ETZERO = E(4) 

CFRAC = S0L(4) 

+++++ 


WRITE OUT FINAL FIT IF REQUESTED 

C+++++ 

IF(IPRFLG(3) .EQ.O) GO TO 250 
SLOPE 1 = SLOPEX(I) 

SL0PE2 = SL0PEX(2) 

SL0PE3 = SL0PEX(3) 

XOFFI = XOFFX(I) 

X0FF2 = X0FFX(2) 

X0FF3 = X0FFX(3) 

TAU1 = -I./SLOPEL 
TAUI = TAU I*TSF 
TAU2 = -I./SL0PE2 
TAU 2 = TAU2*TSF 
TAU 3 = -1./SL0PE3 
TAU 3 = TAU3*TSF 
ETAU1 = ESLL/SL0PE1**2 
ETAUI = ETAUI*TSF 
ETAU2 = ESL2/SLOPE2**2 
ETAU2 = ETAU2*TSF 
ETAU3 = ESL3/SLOPE3**2 
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ETAU3 = ETAU3*TSF 
AREA1 = TAU1*EXP(X0FF1) 

AREA2 = TAU2*EXP(XOFF2 ) 

AREA3 = TAU3*EXP(XOFF3) 

EAREA1 = ETAU1*EXP(EX0L) + EX01*TAU1 

EAREA2 = ETAU2*EXP(EX02) + EX02*TAU2 

EAREA3 = ETAU3*EXP(EX03) + EX03*TAU3 

AREAT = AREA1 + AREA2 + AREA3 
EAREAT = EAREA1 + EAREA2 + EAREA3 
XINT1 = ARE A1 /AREAT 
XINT2 = AREA2 /AREAT 
XINT3 = ARE A3 /AREAT 

E I NT 1 = EAREA1 /AREAT + EAREAT*AREA1 /AREAT**2 
EINT2 = E ARE A2 /AREAT + EAREAT*AREA2 /ARE/ T**2 
EINT3 = EAREA3 /AREAT + EAREAT*AREA3 /ARE/ T**2 
WRITE(6 ,14) SLOPE 1 ,ESLL,SLOPE2 , ESL2.SL01 E3,ESL3, 

1 TAU1 , ETAU1 , TAU2 ,ETAU2 , TAU3 ,ETAU3 , 

2 XOFF1 , EXO L , XOFF2 , EX02 , XOFF3 , EX03 , 

3 AREA1 , EAREA1 , AREA2 , EAREA2 , AREA3 , E ARE A3 . 

4 XINT1 , EINT1 ,XINT2 ,EINT2 ,XINT3 ,EINT3 , 

5 TZERO,ETZERO,STD, ITER1 ,ITER2 ,ICT,IERR 
DELTO = CFRAC - CFRACO 
IF(ICT.NE.O) DELTO = 4.*DELT0 

CFRAC = CFRACO + DELTO 

IF(ABS(CFRAC - CFRACO) .LE. O.OOl) ICT = ICT + 1 
IF (ABS( CFRAC - CFRACO) .GT. 0.001 
1 .OR.ICT.LT. 2) GO TO 390 
GO TO 250 
400 CONTINUE 



C 

C SUBROUTINE CALPLT CLOSES THE PLOT VECTOR FILE. 

C THIS SUBROUTINE IS DOCUMENTED IN "LAIGLEY GRAPHICS 

C SYSTEM", CENTRAL SCIENTIFIC COMPUTIN J COMPLEX 

C DOCUMENT G-3 

C 



CALL CALPLT(0. ,0. ,999) 

STOP 

END 

SUBROUTINE DATPLT(X, Y , ICODE , SLOPE ,XOFF ,N "LA,NPOA,NFIT) 

c 

C SUBROUTINE DATPLT GENERATES GRAPHIC OUTPUT OF THE 

C PROGRAM RESULTS. SEVERAL SUBROUTINES USED BY 

DATPLT ARE DOCUMENTED IN "LANGLEY GRAPHICS 
SYSTEM", CENTRAL SCIENTIFIC COMPUTING COMPLEX 
C DOCUMENT G-3. 

C 

COMMON/PRDAT/DATAX( 800) ,DATAY(800) ,DATAX"(800) ,DATAYL(800) 
COMMON /PEEK/ IPEEK( 2 ) 

C0MM0N/0TDAT/SL0PEX(4) ,X0FFX(4) ,ESL0PE(4 ,EX0FF(4) , 

1 NP1ST(4) ,NPLST(4) ,TAU(4) ,ETAU(4) ,AREA(4 . ,EAREA(4) 

COMMON / 1 NDAT/ICNT( 800) ,NPTS ,NPTSS , BKG, NDJ V , IMAXC ,TSF ,IDCHAN, 
1 IPRFLG(IO) , IPLFLG( 10) ,NP0(4) ,NPL(4) , LABEL, CFRAC 
DIMENSION JMESS(4) ,KMESS(3) 
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o n o o o o 


COMMON /FITDAT/XX( 800) ,YY(800) 

DIMENSION X(800) , Y(800) 

DIMENSION MESS(2) 

DIMENSION MESSI(4) 

DATA IENT/0/ 

DATA JMESS/10HFIRST FIT .10HSECOND FIT , 10HTHIRD FIT , 
1 10HF0URTH FIT/ 

DATA KMESS/IOHOR FIRST ,10H0R SECOND , 10H0R THIRD / 
NPST = NPTS 

IF(IPEEK(1).EQ.0) GO TO 20 
YMAX = 0. 

DO 10 1=1, NPST 

IF( Y( I) . LT .YMAX) GO TO 10 

YMAX = Y(I) 

IMAX = I 
10 CONTINUE 

NPST = IPEEK(2) - IPEEK(l) + 1 
20 CONTINUE 

IF ( I CODE . EQ. 3 .OR . ICODE . EQ. 5) NPST = NPTSS 
I ENT = IENT + 1 


SUBROUTINE NFRAME IS A PART OF THE LANGLEY GRAPHICS 
SYSTEM AND EXECUTES A FRAME ADVANCE 



C — 


IF ( IENT .GT . 1) CALL NFRAME 


SUBROUTINE CALPLT IS A PART OF THE LANGLEY GRAPHICS 
SYSTEM AND MOVES THE PEN TO A NEW LOCATION WITH PEN 
UP OR DOWN. CALPLT CAN ALSO ESTABLISH A NEW PLOT 
ORIGIN 


CALL CALPLT ( 1 . 5 , 1 . 5 , -3) 
MESS(l) = 10HLN (COUNTS 
MESS (2) = 10H) 

IF( ICODE . EQ . 1 ) MESSl(l) = 
IF( ICODE . EQ . 1 ) MESS1(2) = 
IF ( ICODE .EQ .2 ) MESSl(l) = 
IF(IC0DE.EQ.2) MESS1(2) = 
IF(IC0DE.EQ.2) MESS1(3) = 
IF(IC0D£.EQ.3) MESSl(l) = 
IF(IC0DE.EQ.3) MESS1(2) = 
IF(ICODE . EQ . 3) MESS1(3) = 
IF(IC0DE.EQ.3) MESS1(4) = 
IF ( ICODE . EQ . 4 . OR . ICODE . EQ 
IF( ICODE . EQ . 4 .OR . ICODE . EQ 
IF(ICOD£ .EQ. 5) MESSl(l) = 
INDX1 = 1 

IF(IPEEK(1) .NE.O) INDX1 = 
INDX2 = NPST 

IF(IPEEK.(1) .NE.O) INDX2 = 
YSV = DATAYL(INDXl) 


10HRAW SPECTR 
10HUM 

10HSPECTRUM M 
10HINUS BACK.G 
10HR0UND 
10HSPECTRUM C 
10HORRECTED F 
KMESS(NFIT) 

10HFIT 

8) MESSl(l) = 10HFITTED SPE 
8) MESS1( 2) = 10HCTRUM 
JMESS(NFIT) 

IMAX + IPEEK(l) 

IMAX + IPEEK(2) 
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DATAYL(INDXl) = 0. 



C 

C SUBROUTINE ASCALE IS A PART OF THE LANGLEY GRAPHICS 

C SYSTEM AND COMPUTES A SCALING FACTOR FOR AN ARRAY 

C OF DATA. ASCALE ALSO DETERMINES THE DATA MINIMUM 

C 



IF(IENT.EQ.L) CALL ASCALE ( DATAYL ( INDX1 ) , 3. ,NPST , 1 , 10. ) 

DATAYL(INDXL) = YSV 

IF(IENT.EQ.I) YMIN = DATAYL(INDX2+1) 

IF(IENT.EQ.l) YSF = DATAYL ( I NDX2+2) 

Y(INDX2+1) = YMIN 

Y(INDX2+2) = YSF 

ISF = ALOGIO(X(NPST ) / 10 . ) 

SF = L0.**ISF 
VP = X(NPST)/SF 
IVP = VP 
IVP = IVP + I 
VMAX = IVP*SF 
VMAX = VMAX/ 10. 

X(NPST+1) = 0. 

X(NPST+2) = VMAX 
MESS2 = 7HCHANNEL 



C 

C SUBROUTINE NOTATE IS A PART OF THE LANGLEY GRAPHICS 

C SYSTEM AND DRAWS ALPHANUMERIC INFORMATION 

C 

C 

CALL NOTATE (4. 5, 7. 1,0. 2, LABEL , 0 . ,7) 

IF(ICOD£.EQ. 5) GO TO 30 

IF (ICODE.EQ. 4. OR. ICODE.EQ. 8) GO TO 50 



C 

C SUBROUTINE AXES IS A PART OF THE LAf GLEY GRAPHICS 

C SYSTEM AND DRAWS A LINE, ANNOTATES THE VALUE OF A 

C VARIABLE AT SPECIFIED INTERVALS WITT TIC MARKS 

C AND PROVIDES AN AXIS IDENTIFICATION LABEL 

C 



CALL AXES (0 . ,0. ,0. ,10. ,X(NPST+1) ,X(NPSTh 2) , 1. ,0. , 

1 MESS2,0.2,-7) 

CALL AXES (0 . ,0. ,90. ,8. ,Y(INDX2+1) ,Y(INDJ 2+2) ,1. ,0. , 

1 MESS(l) ,0.2,11) 

IF(ICODE.EQ.l) CALL NOTATE ( 4 . 5 , 7 . 5 ,0 . 2 ,M£SS1 ,0 . , 12) 

IF(IC0DE.EQ.2) CALL NOTATE (4 . 5 , 7 . 5 ,0. 2 ,MESS1 ,0 . , 25) 

IF( ICODE . EQ . 3 ) CALL NOTATE (4 . 5, 7. 5, 0.2, MESS 1,0. ,33) 

IF( ICODE . EQ . 4 .OR . ICODE . EQ . 8 ) 

1CALL N0TATE(4 . 5,7 .5,0.2,MESS1,0. ,15) 

IF(IC0DE . EQ. 8) CALL NOTATE (3 .5,7 .5,0.2, >HF INAL,0. ,5) 



C 

C SUBROUTINE LINPLT IS A PART OF THE .ANGLEY GRAPHICS 

C SYSTEM AND DRAWS A LINE BETWEEN AND 'OR DRAWS A NASA 
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onor. non non 


STANDARD SYMBOL AT EACH SUCCESSIVE DATA POINT IN AN 
ARRAY 


C 


CALL LINPLT(X,Y(INDX1) , NPST , 1 ,-l , 22 , 1 , 0) 

30 CONTINUE 

IF(IC0DE.LE.4) GO TO 60 
NP = NPLA + NPOA 
DO 40 1=1, NP 
XX(I) = FLOAT (I - 1) 

YY( I ) = SL0PE*XX( I ) + XOFF 
40 CONTINUE 

XX(NP+1) = 0. 

ISF = AL0G10(XX(NP) / 10 . ) 

SF = 10 . **ISF 
VP = XX(NP) /SF 
IVP = VP 
IVP = IVP + 1 
VMAX = IVP*SF 
VMAX = VMAX/ 10 . 

XX(NP+2) = VMAX 
YY (NP+1 ) = YMIN 
YY(NP+2) = YSF 
XSV1 = X(NP+1) 

XSV2 = X(NP+2) 

YSV1 = Y(NP+1) 

YSV2 = Y(NP+2) 

X(NP+1) = XX(NP+1) 

X(NP+2) = XX(NP+2) 

Y(NP+1) = YY(NP+1) 

Y(NP+2) = YY(NP+2) 

CALL AXES(0. ,0. ,0. ,10. ,X(NP+1) ,X(NP+2) ,1. ,0. , 
1 MESS2 ,0 . 2,-7) 

CALL AXES(0 . ,0. ,90. ,8. ,Y(NP+1) ,Y(NP+2) ,1. ,0. , 
1 MESS(l) ,0.2,11) 

CALL NOTATE (4 .5,7 . 5,0.2,MESS1 ,0. ,10) 

CALL LINPLT(X,Y,NP,1,-1,22,1,0) 

CALL L I N PLT ( XX ,YY,NP, 1,0, 0,0,0) 

XXP = XX(NPOA) /XX(NP+2) 

YYP = YY(NPOA) /YY(NP+2) 


SUBROUTINE PNTPLT IS A PART OF THE LANGLEY GRAPHICS 
SYSTEM AND DRAWS A NASA STANDARD SYMBOL CENTERED ON 
A GIVEN COORDINATE 


CALL PNTPLT (XXP, YYP, 3, 2) 
XXP = XX(NP) /XX(NP+2) 

YYP = YY(NP) /YY(NP+2 ) 
CALL PNTPLT (XXP, YYP, 3, 2) 
X(NP+1) = XSV1 
X(NP+2) = XSV2 
Y(NP+1) = YSV 1 
Y(NP+2) = YSV2 
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GO TO 60 
50 CONTINUE 

Y (NPST+1 ) = YMIN 
Y(NPST+2) = YSF 

CALL AXES(0.,0.,0.,10.,X(NPST+1),X(NPS '+2),1.,0., 

X CALL AXES(0. ,0. ,90. ,8. ,Y(NPST+1) ,Y(NPS C+2) , 1 • ,0 . , 

\f(IC0De’.EQ.8) CALL N0TATE(3 . 5 , 7 . 5 , 0 . 2 , 5HFINAL.0 . , 5) 
CALL NOTATE (4 . 5,7. 5 ,0. 2 ,MESS1 ,0 , 15) 

CALL LINPLT(X,Y,NPST, 1 ,-1,22,1,0) 

YY (NPST+1) = YMIN 

YY(NPST+2) = YSF 

XX( NPST+1) = X(NPST+1) 

XX(NPST+ 2 ) = X(NPST+2) 

CALL LINPLT(XX( 1) , YY( 1) ,NPST , 1 ,0 ,0 ,0 ,( ) 

60 CONTINUE 
RETURN 
END 

SUBROUTINE DATRD 


C 

c 

c 

c 


SUBROUTINE DATRD READS THE EXPERIMENTAL DATA AND 
THE OPERATION INSTRUCTIONS 


COMMON/ INDAT/ICNT (800) ,NPTS ,NPTSS , BKG , NDIV, IMAXC | ,TSF , IDCHAN , 
1 IPRFLG(IO) .IPLFLG(IO) ,NP0(4) ,NPL(4) , .ABEL.CFRAC 
COMMON /MXNIT/ ITERMX 

COMMON /PEER/ IPEER( 2 ) 

COMMON/GAUSDAT/AMP , SIGINV, IGAUSS ,GAUS^(91) IPLFLG 

NAMELIST/ INDAT/BKG, NDIV, TSF , IDCHAN, NPO ,NPL , IPRFLG , IPLF LG , 

1 AMP , SIGINV, IGAUSS , I PEEK, ITERMX, CFRAu 
L FORMAT ( 14, A7) 

REWIND 1 

READ ( 1 , 1) NPTS, LABEL 
READ(1,*) (ICNT(I) ,1=1, NPTS) 

BKG = 0. 


NDIV = 51 

N1 = NPTS - NDIV + 1 
DO 10 I=N 1 , NPTS 
BKG = BKG + FLOAT( ICNT( I) ) 
10 CONTINUE 

BKG = BKG/ FLOAT (NDIV) 

TSF = 61.57 
IDCHAN = -2 
NPO(l) = 120 
NPL(l) = 30 
NP0(2) = 38 
NPL( 2) = 20 
NP0(3) = 11 
NPL(3) = 9 

NPO(l) = NPO(l) - IDCHAN 
NPU(2) = NP0(2) - IDCHAN 
NP0(3) = NP0(3) - IDCHAN 
NP0(4) = 0 
NPL(4) = 0 
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DO 20 1=1,10 
IPRFLG(X) = 0 
IPLFLG(I) = 1 
20 CONTINUE 

IPRFLG(3) = 1 
IPRFLG(6 ) = 1 
IPLFLG(2) = 0 
AMP = 0.22387 
SIGINV = 0.13745 
IGAUSS = 1 
ITERMX = 30 
CFRAC = -0.5 
IPEEK(l) = 0 
IPEEK(2) = 0 
READ( 5 .INDAT , END=30) 

30 IF(E0F ( 5) . NE . 0) GO TO 40 
40 CONTINUE 

IF( IGAUSS . EQ. 0) AMP = 0. 

IF( IGAUSS . NE . 2) GO TO 50 
READ ( 1 , * ) (GAUSS(KK) , KK.=38 ,54) 

50 CONTINUE 

WRITE(6, INDAT) 

RETURN 

END 

SUBROUTINE LSQ2(N , X , Y , SLOPE ,XOFF , ITER , I ERR, STD, El ,£2) 

c 

C SUBROUTINE LSQ2 OBTAINS THE LEAST-SQUARE FIT OF 

C A STRAIGHT LINE TO A PORTION OF THE EXPERIMENTAL 

C DATA IN THE LOG DOMAIN 

C 

DIMENSION X(30) ,Y(30) ,ARAY(2) ,R(100) ,B(2,3) ,C(2) 
DIMENSION KARY( 7 ) ,ERROR( 2) 

DATA ERROR/2*l.E-5/ 

TERR = 0 
XOFF = 0. 

SLOPE =20. 

ITER = 0 
10 CONTINUE 
DO 20 1=1,2 
DO 20 J=1 , 3 
B(I , J) = 0. 

20 CONTINUE 

ITER = ITER + 1 
DO 50 1=1, N 
ARAY(l) = X( I ) 

ARAY(2) = 1. 

FX = SLOPE*X( I ) + XOFF 
R(I) = Y( I ) - FX 
ISUM = 1 
DO 40 K=1 ,2 

B(K, ISUM) = B ( K. , I S UM ) + ARAY(ISUM)**2 
B(K,3) = B(K, 3) + ARAY(K)*R(I) 

IF(ISUM.EQ.2) GO TO 40 
JM1 = ISUM 
ISUM = ISUM + 1 
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DO 30 J=ISUM,2 

B(K,J) = B(K, J) + ARAY(JM1)*ARAY(J) 

30 CONTINUE 
40 CONTINUE 
30 CONTINUE 
DO 70 1=1,2 
IM1 = I - 1 
IF(I.EQ.l) GO TO 70 
DO 60 J=1 , IM1 
B(I,J) = B( J , I ) 

60 CONTINUE 
70 CONTINUE 

KARY(l) = 10 
KARY ( 2 ) = 2 
KARY ( 3 ) = 3 
KARY (4) = 0 
KARY ( 3 ) = 2 
KARY(6) = 0 
KARY ( 7 ) = 0 

c 

C 

C SUBROUTINE MATOPS OBTAINS THE INVERSE OF A MATRIX AND THE 

C SOLUTION OF A SET OF LINEAR EQUATIONS. THIS SUBROUTINE IS 

C DOCUMENTED IN "MATHEMATICAL AND STATISTICAL SOFTWARE AT 

C LANGLEY", CENTRAL SCIENTIFIC COMPUTING COMPLEX DOCUMENT 

C N2-3A 

C 

c 

CALL MATOPS (KARY, B,DET, DUMMY) 

DO 80 1=1,2 
C(I ) = B(I ,3) 

80 CONTINUE 

IF (SLOPE . EQ . 0 . ) GO TO 110 
TEST = ABS(C(1)/SL0PE) 

IF(TEST.GT.ERR0R(1) ) GO TO 150 
90 IF (XOFF . EQ . 0 . ) GO TO 120 
TEST = ABS(C(2) /XOFF) 

IF(TEST.GT.ERR0R(2) ) GO TO 150 
100 CONTINUE 
GO TO 130 

110 IF(ABS(C(1)).GT.ERR0R(1)) GO TO 150 
GO TO 90 

120 IF(ABS(C(2)).GT.ERROR(2)) GO TO 150 
GO TO 100 
130 CONTINUE 

SLOPE = SLOPE + C(l) 

XOFF = XOFF + C(2) 

STD = 0. 

DO 140 I J=1 , N 
STD = STD + R(IJ)**2 
140 CONTINUE 

STD = STD/ FLOAT (N- 2) 

El = B(1,1)*STD 
E2 = B(2 ,2)*STD 
El = SQRT(El) 
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E2 = SQRT(E2) 

STD = FLOAT (N-2)*STD 
STD = SQRT(STD/FLOAT(N) ) 

RETURN 

150 IF( ITER.GT . 30) GO TO 160 
SLOPE = SLOPE + C ( L ) 

XOFF = XOFF + C(2) 

GO TO 10 
160 IERR = 1 
GO TO 130 
END 

SUBROUTINE LSQ6(NPTS ,X ,Y , SOL , E , ITER , IERR , STD , IMAXC , 

1 NPO , NPL , CFRAC , ICT) 

C 

C SUBROUTINE LSQ6 LOCATES THE AMPLITUDES AND LIFETIMES 

C FOR FIXED ZERO TIME THAT MINIMIZE THE STANDARD 

C DEVIATION BETWEEN THE ANALYTIC APPROXIMATION AND THE 

C EXPERIMENTAL DATA 

C 

DIMENSION NPL(4) ,NP0(4) 

DIMENSION FXC(800) ,D£R(6,8Q0) 

DIMENSION X(800) , Y(800) 

COMMON /MXN IT/ ITERMX 

COMMON /GAUS DAT /AMP , S IGINV , IGAUSS ,GAUSS(91) 

DIMENSION E(6) 

DIMENSION R(800) 

DIMENSION SOL( 6) 

DIMENSION ARAY(25) , B( 6 , 7 ) ,C( 25) 

DIMENSION KARY(7) , ERROR/ 25) 

DATA ERROR/25*i.E-5/ 

DATA EPS/l.E-5/ 

DATA NENTRY/O/ 

ITER = 0 
IERR = 99 
E(l) = 1. 

E(2) = 1. 

E(3) = 1. 

E(4) = 1. 

E(5) = 1. 

E(6) = 1. 

NENTRY = NENTRY + 1 
IF(NENTRY.EQ. 1) RETURN 
IERR = 0 
IPT = IMAXC + 3 
IF(ICT.EQ.l) IPT = IMAXC - 6 
ITER = 0 
10 CONTINUE 
DO 20 1=1,6 
DO 20 J=1 , 7 
B( I , J) = 0. 

20 CONTINUE 

ITER = ITER + 1 
NVAR = 6 
NXO = IMAXC 
XO = FLOAT( IMAXC) 
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NPTSP = NPTS + 20 

DO 30 JJ=1, NPTSP 

IF( J J. LT .NXO-16) GO TO 30 

ARGX = FLOAT ( JJ) - (XO + Cf RAC) 

JF(JJ.EQ.NXO) ARGX = 0.5*(0.3 - CFRAC) 

IF(ARGX.LT .0. ) ARGX = 0. 

ARG = S0L(2)*ARGX 
Cl = 0. 

IF(ARG.LT.670. ) CL = EXP(-ARG) 

DER( 1 , JJ) = Cl 

DER(2,JJ) = -S0L(1)*ARGX*DER(1,JJ) 

ARG = S0L(4)*ARGX 
C2 = 0. 

1F(ARG.LT .670. ) C2 = EXP(-ARG) 

DER(3 , JJ) = C2 

DER(4 , JJ) = -S0L(3)*ARGX*DER(3,JJ) 

ARG = S0L(6)*ARGX 
C3 = 0. 

IF(ARG.LT .670. ) C3 = EXP(-ARG) 

DER(5,JJ) = C3 

DER(6,JJ) = -S0L(5) *ARGX*DER( 5 ,JJ) 

FXC(JJ) = SOL( 1)*DER(1,JJ) + SOL(3)*DER(3 , JJ) + S0L(5)*DER(5 ,JJ) 
IF(JJ.EQ.NXO) FXC(JJ) = (0.5 - CFRAC) *FXC( JJ) 

IF(JJ.LT.NXO) FXC(JJ) = 0. 

IF(JJ.LT.NXO) DER(1,JJ) = 0. 

IF(JJ.LT.NXO) DER( 3 , JJ ) = 0. 

IF(JJ.LT.NXO) DER(5 , J J) = 0. 

30 CONTINUE 

DO 80 1=1, NPTS 
R(I) = 0. 

FX = 0. 

DO 40 JJ=1 , 6 
ARAY(JJ) = 0. 

40 CONTINUE 

IF(I.LT.IPT) GO TO 80 
NST = 1-8 
NEND =1+8 
KK = 37 

DO 50 JJ=NST ,NEND 
KK = KK + 1 

ARAY(l) = ARAY(l) + DER( 1 , JJ)*GAUSS(KK) 

ARAY(2) = ARAY ( 2 ) + DER(2 , JJ)*GAUSS(KK) 

ARAY ( 3 ) = ARAY ( 3 ) + DER(3 , JJ) *GAUSS(KK) 

ARAY ( 4 ) = ARAY ( 4 ) + DER(4 , JJ)*GAUSS(KK) 

ARAY ( 5 ) = ARAY(5) + D£R(5 , JJ)*GAUSS(KK) 

ARAY ( 6 ) = ARAY(6) + DER(6 , JJ)*GAUSS(KK) 

FX = FX + FXC( JJ) *GAUSS(KK) 

50 CONTINUE 

R( I ) = Y( I) - FX 
ISUM = 1 
DO 70 K=1 ,NVAR 

B(K, ISUM) = B(K, ISUM) + ARAY( ISUM) **2 
B(K,NVAR+1) = B(K,NVAR+1) + ARAY(K)*R( I) 

IF( ISUM. EQ.N VAR) GO TO 70 
JM1 = ISUM 


25 



ISUM = I SLIM + 1 
DO 60 J=ISUM,NVAR 

B(K,J) = B(K,J) + ARAY ( JM 1 ) *ARAY ( J ) 

60 CONTINUE 
70 CONTINUE 
80 CONTINUE 

DO 100 1=1, N VAR 
IM1 = I - 1 
IF(I.EQ.l) GO TO 100 
DO 90 J= 1 , IM1 
B(I,J) = B(J,I) 

90 CONTINUE 
100 CONTINUE 

KARY(l) = 10 
KARY (2) = 6 
KARY ( 3 ) = 7 
KARY ( 4 ) = 0 
KARY(5) = NVAR 
KARY ( 6 ) = 0 
KARY(7) = 0 

C 

C 

C SUBROUTINE MATOPS OBTAINS THE INVERSE OF A MATRIX AND THE 

C SOLUTION OF A SET OF LINEAR EQUATIONS. THIS SUBROUTINE IS 

C DOCUMENTED IN "MATHEMATICAL AND STATISTICAL SOFTWARE AT 

C LANGLEY", CENTRAL SCIENTIFIC COMPUTING COMPLEX DOCUMENT 

C N2-3A 

C 

c 

CALL MATOPS (KARY, B.DET, DUMMY) 

DO 110 1=1, NVAR 
C(I) = B(I .NVAR+l) 

110 CONTINUE 

DO 120 K=1 , NVAR 
TST = ABS(C(K) ) 

IF(ABS(SOL(K) ) .GT .EPS) TST = ABS(C(K) /SOL(K) ) 
IF(TST.GT.ERROR(K) ) GO TO 170 
120 CONTINUE 

DO 130 IJ=1 ,NVAR 
SOL(IJ) = SOL(IJ) + C(IJ) 

130 CONTINUE 
140 CONTINUE 
STD = 0. 

DO 150 IJ=1 ,NPTS 
IF(IJ.LT.IPT) GO TO 150 
STD = STD + R(IJ)**2 
150 CONTINUE 

STD = STD/FL0AT(NPTS-NX0-2) 

DO 160 IJ=1 ,NVAR 

E( I J) = SQRT(ABS(B(IJ, IJ))*STD) 

160 CONTINUE 

STD = FL0AT(NPTS-NX0-2)*STD 
STD = SQRT(SID/FLOAT(NPTS-NX0) ) 

RETURN 

170 IF ( ITER.GT . ITERMX) GO TO 200 
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RATMAX = 0. 

DO 180 IJ=1 ,NVAR 

RAT = ABS(C(IJ)/SOL(IJ)) 

IF ( RAT. GT. RATMAX) RATMAX = RAT 
180 CONTINUE 
FRAC = 1. 

I F ( RATMAX . GT . 0 . 4 ) FRAC = 0.4/ RATMAX 

DO 190 IJ=1 ,NVAR 

SOL(IJ) = SOL(IJ) + FRAC*C( IJ) 

190 CONTINUE 
GO TO 10 
200 IERR = 1 
GO TO 140 
END 

SUBROUTINE LSQl 

SUBROUTINE LSQl LOCATES THE VALUE OF ZERO TIME FOR 
FIXED AMPLITUDES AND LIFETIMES WHICH MINIMIZES THE 
STANDARD DEVIATION BETWEEN THE ANALYTIC APPROXIMATION 
AND THE EXPERIMENTAL DATA 

COMMON /MXN IT/ ITERMX 

COMMON /FUNC0M/R( 800) ,X0,NX0 

COMMON/GAUSDAT/AMP , SIGINV , I GAUSS ,GAUSS(9 t ) 

DIMENSION E( 7 ) 

DIMENSION S0L(7 ) 

COMMON/NEWCOM/SOL.E, ITER, IERR, STD 

COMMON /PRDAT/DATAX(800) ,DATAY(800) ,DATAX 1(800) ,DATAYL(800) 
COMMON /CLCDAT/X( 800) ,Y(800) ,YL(800) 

COMMON /INDAT / ICNT(800) , NPTS , NPTSS , BK.G, ND IV , IMAXC ,TSF , IDCHAN , 
1 IPRFLG(IO) ,IPLFLG(10) ,NP0(4) ,NPL(4) , LABEL, CFRAC 
DATA EPS/1. E-5/ 

DATA NENTRY/O/ 

NENTRY = NENTRY + 1 

IERR = 0 

ITER = 0 

NXO = IMAXC 

XO = FLO AT (IMAXC) 

CALL FUNX(STDO.SOL) 

DELX =0.05 

IF (NENTRY. GT.l) DELX = 0.01 
SOL(4) = SOL(4) + DELX 
CALL FUNX(STD , SOL) 

IF(STD.GE.STDO) DELX = -DELX 
STDO = STD 

SOL(4) = SOL(4) + DELX 
10 CONTINUE 

ITER = ITER + 1 
CALL FUNX(STD , SOL) 

IF(STD.GE.STDO) GO TO 20 
STDO = STD 

S0L(4) = S0L(4) + DELX 
GO TO 10 

20 TST = ABS (STDO - STD) /STDO 
IF(TST.LT.EPS) GO TO 30 


27 



S0L(4) = S0L(4) - 2 . 0*DELX 
DELX = 0 . 5*D£LX 
STDO = 1.E99 
GO TO 10 
30 CONTINUE 
STD = 0. 

DO 40 I J=1 , NETS 
IF(IJ.LT.NXO) GO TO 40 
STD = STD + R(IJ)**2 
40 CONTINUE 

STD = STD/ FLOAT (NPTS-NXO-2) 

E(4) = ABS(DELX) 

RETURN 

END 

SUBROUTINE FUNX(STD, SOL) 

C 

C SUBROUTINE FUNX CALCULATES THE STANDARD DEVIATION 

C BETWEEN THE ANALYTIC APPROXIMATION AND THE EXPERIMENTAL 

C DATA. THE THREE AMPLITUDES ARE IN SOL(I), SOL(2), 

C AND SOL(3) . THE ZERO TIME IS IN S0L(4). THE THREE 

C LAMBDAS (1/TAU) ARE IN SOL(5), S0L(6), AND SOL(7). 

C 

COMMON/FUNCOM/R(800) ,X0,NX0 

COMMON /PRDAT/DATAX( 800) ,DATAY(800) ,DATAXT(800) ,DATAYL(800) 
DIMENSION S0L(7 ) 

COMMON/ GAUSDAT/AMP , SIGINV , I GAUSS ,GAUSS(9 1) 

COMMON / INDAT / ICNT ( 800 ) ,NPfS ,NPTSS , BK.G,NDIV,IMAXC,TSF ,IDCHAN , 
1 IPRFLG(IO) ,IPLFLG(10) ,NP0(4) ,NPL(4) .LABEL, CFRAC 
DIMENSION CONT(800) 

NPTSP = NPTS + 20 

DO 10 JJ=1, NPTSP 

IF( JJ . LT .NX0-20) GO TO 10 

ARG = SOL ( 5 ) * ( FLOAT ( J J ) - (XO + S0L(4))) 

IF(JJ.EQ.NXO) ARG = SOL( 5) *0 . 5* (0 . 5 - S0L(4)) 

Cl = 0. 

IF (ARG. LT . 670 . ) Cl = EXP(-ARG) 

ARG = S0L(6)*(FL0AT( JJ) - (XO + S0L(4))) 

IF(JJ.EQ.NXO) ARG = S0L( 6 ) *0 . 5* (0 . 5 - S0L(4)) 

C2 = 0. 

IF(ARG.LT .670. ) C2 = EXP(-ARG) 

ARG = S0L( 7 )* (FLOAT ( J J ) - (XO + S0L(4))) 

IF(JJ.EQ.NXO) ARG = SOL( 7 ) *0 . 3* (0 . 5 - S0L(4)) 

C3 = 0. 

IF(ARG.LT.670. ) C3 = EXP(-ARG) 

CONT ( JJ) = S0L(1)*C1 + SOL(2)*C2 + SOL(3)*C3 
IF(JJ.EQ.NXO) CONT(JJ) = (0.5 - S0L(4) )*CONT( JJ) 
IF(JJ.LT.NXO) CONT(JJ) =0. 

10 CONTINUE 

DO 30 1=1, NPTS 
FX = 0. 

IF(I.LT.NX0-8) GO TO 30 
NST = 1-8 
NEND =1+8 
KK = 37 

DO 20 JJ=NST,NEND 
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KK. = KK. + 1 

FX = FX + CONT(JJ) *GAU S S ( KK. ) 
20 CONTINUE 

R( I ) = DATAY(I) - FX 
30 CONTINUE 
STD = 0. 

DO 40 IJ-l.NPTS 
I F( I J . LT . NXO-8) GO TO 40 
STD = STD + R(IJ)**2 
40 CONTINUE 
RETURN 
END 
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TOTAL 

COUNTS = 

1872471 



NUMBER 

OF CHANNELS 

IN DATA 

FILE = 

650 

BACKGROUND AVG. = 

1.2902 OE +01 

NO I V * 


NUMBER OF COMPONENTS TO FIT * 3 

Figure Al. Input spectrum. 
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RAW SPECTRUM 

CHANNEL TIME (PS) 


1. 

6. 15700E+01 

2 • 

1. 23140E+02 

3 • 

1. 84710E+02 

4, 

2 .46280E+02 

5 • 

3 • 07o5 uE +0 2 

6 • 

3 .69420E+02 

7 • 

4. 30990E+02 

8 . 

4 . 9256OE+02 

9 . 

5.54130E+02 

10 . 

6.157C0E+02 

11 • 

6.77270E+02 

12 • 

7. 38840E+02 

13. 

8. 00410E + 02 

1 A. 

8.61980E+02 

15 . 

9.2355OE+02 

16 . 

9.8512CE+02 

17 • 

1.04669E+03 

18. 

1. 10826E + 03 

19 . 

1.1698 3E + 03 

20 . 

1.23140E+03 

21 • 
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22. 

1 • 354 5<tE + 0 3 

23 . 

1. 4 16 1 1 E+ 03 

24 • 

1.4776oE+03 

25. 

1 .53925E+03 

26. 

1 .60082E+03 

27 • 

1.662396+03 

28. 

1. 72 39&E + 03 

29. 
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30 . 

1.84 710E + 03 

31. 
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32. 

1.97024E+03 
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6.00000E+00 
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I. 10000E+01 
9. OOOOOE+OO 


Figure A2. Sample raw spectrum output. 
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1*608 53 E+ 05 


SPECTRUM MAXIMUM * 

CHANNEL NUMBER CORRESPONDING TO SPECTRUM MAXIMUM 

CORRECTED CHANNEL * 103 

Figure A3. Sample peak output. 


FIRST FIT (FROM 137 TO 167) 

SLOPE * -3 *54o9o816E-02 ( 3 . Ilo9<*633£-03 > 

TAU = 1 . 7 35 8 H8 68 E+0 3 ( 1.52540054E+02) 

INTERCEPT = 9.46815652E+00 ( <.*69tt7 55b7E-01 ) 

AREA = 2.2*636328E+07 ( 1 .9 7483 780E +Oo ) 

STANDARD DEVIATION * 1 .42757049E-01 

ITERATIONS * 2 

I ERR = 0 

Figure A4. Sample fit output (initial estimate). 
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FITTED SPECTRUM 


CHANNEL 

TIME (PS) 

• 

• 

• 

97. 

• 

• 

• 

5.97229E+03 

98 . 

6.0338&E+03 

99. 

6.09543E+03 

100 . 
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101 . 
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102 . 
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124. 
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126. 

7.75782E+03 

127 . 
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128 . 

7.8 8 09oE +0 3 

129. 

7.94253E+03 

130. 

• 

• 

8. 00410 E + 03 

• 

• 

• 

NONLINEAR 

• 

FIT STANDARD DEVIATI 

VARIANCE 

* 1.80950085E+00 

Figure A5. Sample fitted 
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Appendix B 
Least-Squares Solution 

The solution of equation (Al) consists of seven 
parameters: A \ , A 2) A3, Ai, A 2 , A3, and t z , where 
A x is the amplitude of the zth component, t z is time 
zero, A z — 1 / , and r x is the lifetime of the ith 
component. The method of least squares has been 
used to obtain the solution. This method is ex- 
tremely sensitive to the initial estimated values of 
the seven parameters; accordingly, a multistep pro- 
cedure is employed to obtain crude initial estimates 
and then refine them. This procedure is described in 
appendix A. 

In general, equation (Al) can be written as 
follows: 

n k ~ B k (tj, Ai, A 2 , A 3 , X 2 , + £fc (Bl) 

where the measurement error in channel k is given 
by Expanding equation (Bl) in a Taylor series 
and dropping higher order terms yields the following 
linear approximation: 

n k ~ ( n k ) 0 + ^ (Ai - Aifi) -f (A 2 - A 2 } 0) 

+ |g(^-A 3 ,„) + ^(A I -A 1 , 0 ) 

+ |g( A2 - A2 , o) 

dFjr , , 

+ d \ J ( A 3 - ^°) + £/t (B 2 ) 

Equation (B2) can be simplified with the following 
definitions: 



Aai = Al - Al'O Aq 2 = A2 - ^ 2,0 Aa 3 = A3 - ^3,0 
A04 — Ai - Ai ? o Aa5 = A2 — A2 ,o Aa6 = A3 — Aa^ 

The simplified equation (B2) takes the form: 

6 

a n k = E b k,i^O‘l + (B 3 ) 

1=1 

or, in matrix notation, 

z k = B k Aa + e k (B 4 ) 


Z k = An*., 

' Aai ‘ 

Ac*2 

Aa = , 

.A a 6 _ 

B k = b k, 2 > b k, 3 > b kfi] 1 

and 

= [ek] 

Thus, for N channels, there are N matrix equations, 
which may be written as 

Z = BAa + e (B5) 

where 



The sum of the squares of the residuals 

e T e= (Z - BAa) T (Z - BAa) (B6) 

is a minimum for that value of Aa for which the 
first variation with respect to Aa of equation (B6) 
vanishes if the second variation with respect to Aa 
is positive. Both these conditions are satisfied for 

A a'=(B T By l B T Z (B 7 ) 

The value of Aa obtained from equation (B 7 ) is 
used to determine new values of A % and A;, and 
this iteration process continues until Aa' — ► 0 . This 
process leads to the best values of A x and A^. 

The derivatives included in equation (B 2 ) are 
readily determined from equation (Al) and are given 
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by the following expressions: 


dFk 

dA\ 


= J2 w 3A e ~ Xl{tj ~ tz) 

3=0 


dFk 
dA 2 


= T, W J,d e ~ X2{tj ~ tz) 


3=0 


dFk 

dA 3 


E 

3=0 


Wi 


k e 


'•^3 


dF\ 

d\ x 


= ~E W 3^ t 3~ t M ie 


— Aj ( tj—tz ) 


3=0 


dF\ 

dX 2 


= -E W 3,k{tj-tz)Me- X * {t ’- tz) 

3=0 


dfk 

dX 3 


E W 3,k( l 3 - ^) A S e 


3=0 


-^3 {tj-tz) 


(B8) 
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Appendix C 

Positron Lifetime Measurements in Epoxy 
Samples 

As part of an ongoing study for elucidating the 
role of metal ions in controlling free volume in poly- 
mers, the following chromium compounds were intro- 
duced in Epon 828 epoxy resin. The metal-complex 
concentration in each case was 1 for every 10 repeat 
units of the host epoxy. 

1. Chromic acetate [Cr(Ac) 3 ] 

2. Chromic perchlorate [Cr(DMSO)e ( 0104 ) 3 ] 

3. Chromous chloride [CrCl 2 ] 


The samples were prepared in the form of 1- in- 
diameter by 0.1-in-thick discs. Positron lifetime spec- 
tra were measured in these samples by using a stan- 
dard fast-fast coincidence measurement system. (See 
refs. 12 to 15.) The experimental lifetime spectra 
are listed in table CI. These spectra have been ana- 
lyzed using the program for analyzing positron life- 
time spectra (PAPLS) and POSITRONFIT program. 
The results of these analyses have been summarized 
in tables I to III in the text. 


38 



Table Cl. Experimiental Lifetime Spectr i Observed in Epon 828 Epoxy 
Sample Containing 0.1 Mole Fraction of Transition Metal Compounds 


(a) Cr(Ac) 3 


85640 

0 

0 

1 

2 

0 

2 

4 

6 

6 

5 

3 

6 

12 

12 

7 

4 

9 

6 

11 

7 

10 

5 

12 

51 

235 

91292 

95598 

35696 

30604 

12580 

11396 

5811 

5439 

3525 

3121 

2375 

2126 

1626 

1470 

1104 

1091 

801 

783 

620 

611 

497 

406 

345 

308 

259 

248 

218 

177 

156 

141 

103 

115 

88 

68 

61 

61 

58 

55 

36 

35 

36 

38 

17 

24 

15 

29 

23 

23 

18 

10 

19 

12 

16 

15 

13 

18 

8 

11 

10 

13 

19 

10 

8 

8 

11 

9 

6 

13 

8 

6 

9 

10 

8 

10 

6 

14 

3 

9 

4 

9 

6 

3 

6 

6 

9 

6 

6 

10 

7 

6 

7 

12 

8 

8 


0 

0 

2 

1 

2 

0 

5 

6 

3 

4 

3 

4 

8 

6 

11 

16 

6 

10 

9 

9 

11 

7 

9 

6 

939 

3467 

90162 

78976 

26904 

23266 

10186 

9117 

4998 

4614 

3131 

2969 

2079 

2027 

1406 

1311 

1042 

927 

761 

736 

555 

531 

413 

424 

327 

303 

237 

214 

187 

184 

138 

127 

116 

93 

86 

71 

73 

53 

47 

46 

39 

52 

36 

31 

28 

35 

22 

28 

21 

17 

21 

20 

15 

11 

14 

13 

20 

13 

12 

7 

11 

10 

5 

7 

13 

11 

11 

9 

9 

9 

11 

6 

8 

10 

11 

13 

5 

10 

10 

5 

6 

8 

9 

2 

10 

8 

11 

9 

5 

9 

13 

13 

8 

8 

11 

7 


0 

0 

1 

1 

5 

4 

5 

2 

5 

3 

9 

7 

7 

10 

5 

11 

10 

5 

16 

9 

5 

8 

14 

9 

11143 

25965 

67947 

57663 

20652 

18178 

8407 

7534 

4314 

3947 

2706 

2619 

1915 

1787 

1325 

1247 

963 

932 

690 

694 

534 

471 

410 

390 

293 

273 

220 

227 

180 

191 

137 

110 

113 

94 

80 

84 

61 

57 

51 

30 

41 

29 

20 

29 

36 

22 

25 

26 

19 

22 

18 

11 

16 

12 

13 

10 

14 

15 

13 

10 

11 

9 

10 

9 

4 

9 

8 

6 

7 

9 

12 

8 

19 

5 

9 

9 

9 

8 

7 

10 

5 

10 

6 

12 

10 

10 

6 

8 

9 

6 

7 

7 

8 

7 

30 

5 


0 

0 

1 

1 

0 

0 

6 

3 

4 

4 

6 

3 

9 

9 

6 

10 

7 

9 

7 

6 

6 

8 

16 

22 

48945 

73446 

49008 

41536 

15907 

14229 

6821 

6286 

3753 

3587 

2528 

2328 

1698 

1573 

1142 

1160 

859 

847 

639 

626 

485 

493 

405 

309 

294 

264 

180 

205 

172 

158 

140 

131 

95 

95 

75 

67 

67 

54 

43 

50 

31 

29 

31 

26 

22 

24 

23 

19 

19 

26 

26 

21 

14 

17 

16 

6 

12 

17 

11 

17 

8 

13 

11 

5 

3 

10 

8 

6 

12 

7 

10 

5 

9 

5 

8 

5 

2 

8 

8 

6 

10 

11 

2 

6 

7 

6 

8 

8 

10 

9 

9 

14 

8 

12 

11 

10 
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Table CL Continued 


(b) Cr(DMSO) 6 (C10 4 ) 3 


85364 

0 

0 

0 

2 

3 

1 

2 

1 

4 

3 

4 

2 

5 

5 

6 

5 

2 

4 

3 

7 

3 

7 

7 

8 

5 

10 

7 

10 

8 

9 

11 

5 

9 

11 

4 

8 

7 

8 

15 

9 

14 

8 

7 

10 

10 

6 

14 

77 

242 

972 

3619 

96094 

102257 

96753 

85425 

38232 

32262 

28756 

24449 

12941 

11659 

10281 

9238 

5642 

5161 

4800 

4389 

3208 

2853 

2807 

2635 

1932 

1954 

1848 

1728 

1377 

1277 

1253 

1186 

955 

958 

890 

897 

683 

711 

656 

639 

554 

479 

501 

456 

410 

380 

395 

362 

315 

295 

269 

250 

202 

204 

181 

179 

175 

145 

163 

164 

114 

123 

107 

115 

107 

110 

95 

89 

68 

75 

75 

67 

52 

58 

49 

58 

50 

57 

45 

37 

33 

31 

33 

36 

34 

38 

27 

31 

19 

25 

14 

23 

19 

18 

20 

34 

26 

27 

16 

26 

16 

14 

16 

8 

11 

9 

15 

10 

10 

11 

15 

14 

14 

11 

8 

14 

10 

12 

10 

13 

9 

10 

12 

14 

17 

5 

13 

11 

10 

7 

5 

10 

6 

5 

7 

8 

6 

8 

10 

9 

7 

9 

13 

14 

10 

9 

10 

3 

6 

10 

10 

8 

6 

5 

5 

13 

14 

10 

4 

2 

5 

8 

11 

15 

10 

3 

13 

11 

11 

10 

9 

7 

9 

3 

10 

8 

12 

8 

4 

2 

7 

7 

6 

7 

12 

4 

9 

4 

11 

10 

6 

9 


0 

0 

0 

1 

0 

1 

1 

0 

5 

0 

3 

2 

8 

6 

5 

3 

6 

5 

6 

5 

1 

8 

4 

4 

5 

8 

13 

15 

7 

9 

9 

6 

7 

9 

9 

6 

4 

12 

11 

11 

6 

10 

7 

12 

7 

10 
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Table I. Summary of Input Parameters 
F = Lifetime of ith component 1 


I{ = Intensity of ith component 


\_t z — Time zero 

j 


Case 1 

Case 2 

T\ = 300 psec; I\ = 70% 

r\ = 300 psec; I\ = 85% ^ 

T 2 = 700 psec; I 2 = 15% 

T 2 = 700 psec; / 2 = 10% 

r 3 = 2100 psec; 7 3 = 15% 

T 3 = 2100 psec; / 3 — 5% 

t z = Channel 50.300 
Background counts = 10.00 



Table II. Summary of Results Obtained for Case 1 

T{ = Lifetime of ith component 
I x — Intensity of ith component 
t z ~ Time zero 


PAPLS results 

POSITRONFIT results [ 

n = 298 ± 1 psec; h = 70.49 ± 0.02% 

T 2 = 672 ± 16 psec; / 2 = 15.14 ± 0.02% 

r 3 = 2080 ± 22 psec; / 3 = 14.37 ± 0.02% 

n = 310 ± 2 psec; /j = 73.80 ± 0.89% 
t 2 = 771 ± 39 psec; I 2 = 12.39 ± 0.67% 
r 3 = 2128 ± 19 psec; / 3 = 13.81 ± 0.33% 

t z = Channel 50.320 
Background counts = 10.27 
Standard deviation = 27.46 

t z — Channel 50.285 
Background counts = 9.89 
Standard deviation — 85.06 


Table III. Summary of Results Obtained for Case 2 

r x = Lifetime of zth component 
I-i = Intensity of zth component 
t z = Time zero 


PAPLS results 

POSITRONFIT results 

Ti = 298 ± 1 psec; /j = 84.44 ± 0.04% 
r 2 = 643 ± 23 psec; / 2 = 10.52 ± 0.03% 
r 3 = 2021 ± 59 psec; / 3 = 5.04 ± 0.02% 

n = 308 ± 1 psec; T t = 87.64 ± 0.64% 
r 2 = 778 ± 42 psec; / 2 = 7.95 ± 0.49% 
r 3 = 2144 ± 41 psec; 7 3 = 4.45 ± 0.22% 

t z — Channel 50.321 
Background counts = 10.27 
Standard deviation = 27.64 

t z = Channel 50.288 
Background counts = 9.87 
Standard deviation = 98.97 




Table IV. Comparison of POSITRONFIT and PAPLS Programs for Lifetime Spectra Analysis 
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Channel 103.11 t z = Channel 103.11 




In(counts) 


16 r 



Channel number (1 channel = 61. 23 psec) 

Figure 1. Computer-generated spectrum with background counts added. 


46 





In(counts) 



Channel number (1 channel 1 61. 23 psec) 

Figure 2. Computer-generated spectrum that shows effect of finite resolution of counting system. 
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In (counts ) 


16 1 — 



Channel number (1 channel = 61. 23 psec) 

(a) Composite spectrum. 

Figure 3. Typical positron lifetime spectrum in an epoxy target. 
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Figure 3. Concluded. 
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16 


14 



Channel number (1 channel * 61. 23 psec) 


(a) Experimental spectrum, minus background, showing third component; t z = time zero. 


Figure 4. Analysis procedure for a positron lifetime spectrum. 
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In (counts) 


16 



Channel number (1 channel - 61. 23 psec) 


(b) Experimental spectrum, minus background and third component, showing second component clearly. 

Figure 4. Continued. 
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In(counts) 


16 


14 


12 



Channel number (1 channel = 61. 23 psec) 

(c) Experimental spectrum showing second component fit. 


Figure 4. Continued. 
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Figure 4. Continued. 
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In(counts) 


16 




Channel number (1 channel = 61. 23 psec) 

(e) Residual spectrum showing first component fit. 
Figure 4. Continued. 
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